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Abstract 

Studies of the strangeness changing hadronic weak interaction have produced a 
number of puzzles that have so far evaded a complete explanation within the Stan- 
dard Model. Their origin may lie either in dynamics peculiar to weak interactions 
involving strange quarks or in more general aspects of the interplay between strong 
and weak interactions. In principle, studies of the strangeness conserving hadronic 
weak interaction using parity violating hadronic and nuclear observables provide 
a complementary window on this question. However, progress in this direction 
has been hampered by the lack of a suitable theoretical framework for interpreting 
hadronic parity violation measurements in a model-independent way. Recent work 
involving effective field theory ideas has led to the formulation of such a framework 
while motivating the development of a number of new hadronic parity violation 
experiments in few-body systems. In this article, we review these recent develop- 
ments and discuss the prospects and opportunities for further experimental and 
theoretical progress. 
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1 Introduction 



Explaining the weak interactions of quarks in terms of the dynamics of the Standard 
Model (SM) has been an area of vigorous research in nuclear and particle physics for 
several decades. Experimentally, the hadronic weak interaction (HWI) is probed by 
observing non-leptonic, flavor changing decays of mesons and baryons and by measur- 
ing observables that conserve flavor but violate the parity symmetry of the strong and 
electromagnetic interactions. Theoretically, the problem has been a particularly chal- 
lenging one, requiring the computation of low-energy weak matrix elements of the HWI 
in strongly interacting systems. Although the structure of the weak quark-quark inter- 
action in the SM has been well established for some time, its manifestation in strongly 
interacting systems remains only partially understood. The stumbling block has been 
the non-perturbative nature of quantum chromodynamics (QCD) at low energies. In 
contending with it, theorists have resorted to a variety of approximation schemes to ob- 
tain physically reasonable estimates of HWI observables. Ultimately, however, arriving 
at definitive, SM predictions requires that one treat the non-perturbative QCD dynamics 
in a rigorous way. 

In the case of the flavor changing decays of mesons, use of effective field the- 
ory (EFT) techniques - chiral perturbation theory (%PT), heavy quark effective the- 
ory (HQET), and recently, soft collinear effective theory (SCET) - have led to enormous 
progress. In each instance, the presence of distinct physical scales at play in the processes 
of interest allows one to carry out a systematic expansion of the effective Lagrangian in 
powers of scale ratios while incorporating the symmetries of QCD into the structure 
of the operators. The operator coefficients that encode the non-perturbative QCD dy- 
namics are obtained from measurement, and the structure of the EFT is then used to 
translate this information into predictions for other observables. Moreover, a meaning- 
ful confrontation of experiment with QCD theory can be made, as computations of the 
operator coefficients can in principle be performed on the lattice. 

The situation involving the HWI of baryons is far less satisfactory, and decades of 
experimental and theoretical work have left us with a number of unresolved puzzles. In 
the case of hyperon non-leptonic decays, for example, one has not yet been able to find a 
simultaneous accounting of both the parity conserving P-wave and parity violating (PV) 
S-wave decay amplitudes. Similarly, the PV asymmetries associated with the radiative 
decays of hyperons are anomalously large. In the limit of degenerate u-, d-, and s- 
quarks, SU(3) flavor symmetry implies that these asymmetries must vanish. Given the 
known mass splitting between the strange and two light flavors, one would expect the 
asymmetries to have magnitudes of order m s /MB ~ 0.15, where Mb ~ 1 GeV is a 
typical hyperon mass. The experimental asymmetries, in contrast, are four-to-five times 
larger in magnitude. Even the well-known AI =1/2 rule that summarizes the observed 
dominance of the 1=1/2 channel over the 7 = 3/2 channel in strangeness changing 
nonleptonic decays remains enigmatic, as no apparent symmetry favors either channel. 
In short, consideration of QCD symmetries and the relevant physical scales does not 
suffice to account for the observed properties of the AS* = 1 HWI. 

While the puzzles surrounding the strangeness changing HWI have been discussed 
extensively elsewhere, the AS = HWI has generally received less attention. Nonethe- 
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less, since we do not know whether the breakdown of QCD symmetry-based expectations 
in the AS = 1 sector results from the presence of a dynamical strange quark or from 
other, yet-to-be-uncovered dynamics, consideration of the AS = HWI - for which the 
strange quark plays a relatively minor role - is no less important. In the following review, 
we focus on this component of the HWI. 

According to the SM, the structure of the low-energy AS = HWI is relatively 
simple: 



where Gp is the Fermi constant and where J^ c and J^ c are the weak charged and 
neutral currents, respectively. The theoretical challenge is to find the appropriate effec- 
tive interaction Hhwi (N, ir, A, . . .) that best describes the hadronic manifestation of 
7^hwT°- Because Jy transforms as a doublet under strong isospin while Jf c contains 
1 = and 1 = 1 components, the current- current products in T^hwi contain terms that 
transform as isoscalars, isovectors, and isotensors. Consequently, T^hwi eS mu st contain 
the most general set of operators having the same isospin properties. In what follows, 
we review the theoretical efforts to determine this effective interaction. 

Experimentally, the AS = HWI can be isolated solely via hadronic and nuclear 
physics processes that violate parity, thereby filtering out the much larger effects of the 
strangeness conserving strong and electromagnetic interactions. Efforts to do so are 
not new. Soon after the 1957 discovery of parity violation in /i-decay and nuclear (3- 
decay, the search was on for evidence of a PV weak nuclear force that would result in 
small, parity violating effects in nuclear observables. That year, Tanner reported the 
first experimental search for a PV nucleon-nucleon (NN) interaction (0). Subsequently, 
Feynman and Gell-Mann (0) predicted that the four fermion interactions responsible 
for leptonic and semi-leptonic weak decays should have a four nucleon partner that is 
similarly first order in Gf- A decade later, Lobashov et al. produced the first definitive 
evidence for the existence of a first order weak NN force in radiative neutron capture on 
181 Ta that was consistent with the Feynman and Gell-Mann hypothesis (0,0). 

The pursuit of this evidence in the Tanner, Lobashov and subsequent experiments 
was challenging, as one expected the magnitude of the PV effects to be 0(1CT 7 ). Along 
the way, it was realized that certain accidents of nuclear structure in many-body nuclei 
could amplify the expected PV effects by several orders of magnitude, and a ~ 10% PV 
effect was, indeed, observed in 139 La (0). The amplification arises from two sources: the 
presence of nearly degenerate opposite parity states that are mixed by the HWI, and the 
interference of an otherwise parity forbidden transition amplitude with a much larger 
parity allowed one. Subsequent experiments then yielded a mix of PV measurements in 
nuclei, where one expected amplification factors of order 10 2 to 10 3 , as well studies of PV 
observables in the scattering of polarized protons and neutrons from hadronic targets. 

Theoretically, however, the use of nuclear systems introduces an additional level of 
complication in the interpretation of experiments, as one must contend with both nu- 
clear structure effects as well as the dynamics of non-perturbative QCD. For over two 
decades now, the conventional framework for carrying out this interpretation has been a 
meson exchange model, popularized by the seminal work of Desplanques, Donoghue, and 
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Holstein (DDH) (6). The model assumes that the PV nucleon-nucleon (NN) interaction 
is dominated by the exchange of the pion and two lightest vector mesons (p and uj), 
and its strength is characterized by seven PV meson-nucleon couplings: h\, hp' 1 ' 2 , h ' 
and h^j 1 , where the superscript indicates the isospin 1 . DDH provided theoretical "rea- 
sonable ranges" and "best values" for the h l M using SU(6) symmetry, constraints from 
non-leptonic hyperon decay data, and the quark model to estimate the experimentally 
unconstrained terms. Despite various attempts to improve upon the original DDH work, 
the results of their analysis still remain as the benchmark, theoretical targets for the PV 
meson-nucleon couplings. 

The experimental results from nuclear and hadronic PV measurements have been 
analyzed using the DDH framework, leading to constraints on combinations of the h l M 
that typically enter PV observables. The results are in general agreement with the DDH 
reasonable ranges, though the ranges themselves are quite broad, and the constraints from 
different experiments are not entirely consistent with each other. A particular quandary 
involves h\\ the 7-decays of 18 F imply that it is consistent with zero, while the analysis 
of the 133 Cs anapole moment differs from zero by several standard deviations ®. More 
to the point, the connection between the PV experiments and SM expectations is far 
from transparent. Indeed, in order to draw this connection using the meson-exchange 
framework and nuclear PV observables, one has to sort through a number of model 
dependent effects involving nuclear structure, hadron structure, and the meson exchange 
model itself. Whether one has a reasonable hope for doing so in a systematic manner is 
debatable at best. 

At the end of the day, the goal of studying the AS = HWI with hadronic and 
nuclear PV is to help determine the degree to which the symmetries of QCD characterize 
the realization of the HWI in strongly interacting systems and, as a corollary, to shed 
light on the long standing puzzles in the AS = 1 sector. To that end, one would ideally 
formulate the problem to make the contact with the underlying SM as transparent as 
possible while avoiding hadronic model and nuclear structure ambiguities. Recently, a 
framework for doing so has been formulated in Reference IQ) using effective field theory 
ideas. That work builds on the extensive developments in the past decade of an EFT 
for the strong NN interaction that has been applied successfully to a variety of few-body 
nuclear phenomena. In the case of the PV NN force, two versions of the EFT are useful, 
depending on the energy scales present in the process under consideration: 

(I) For energies well below the pion mass, the EFT contains only four-nucleon operators 
and five effective parameters, or "low-energy constants" , that characterize the five 
independent low-energy S-wave/P-wave mixing matrix elements: A^' 1 ' 2 , X t , and 
pt- Relative to the leading order parity-conserving four nucleon operators, the PV 
operators are 0(Q), where Q is a small energy scale. In this version of the EFT, 
the pion is considered to be heavy and does not appear as an explicit, dynamical 
degree of freedom. 

(II) At higher energies, the the pion becomes dynamical and three additional constants 

1 In the literature, the isovector, PV irNN coupling is often denoted / w . Here, however, we adopt the 
h\ notation to avoid confusion with the pion decay constant. 
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associated with 7r-exchange effects appear at lowest order: h^, along with a second 
parameter in the EFT potential, k\ a , and a new meson-exchange current opera- 
tor characterized by C n . Moreover, the EFT incorporates the effects of two-pion 
exchange for the first time in a systematic way, leading to predictions for a medium- 
range component of the PV NN interaction. 

The essential differences between the PV EFT and the meson-exchange frameworks - as 
well as their similarities - are summarized in Figure ^ 




(a) (b) (c) (d) 



Figure 1: Comparison of (a) meson- exchange and (b-d) effective field theory (EFT) treat- 
ments of the parity-violating NN interaction. Panels (b), (c), and (d) give illustrative 
contributions to short, medium, and long-range components, respectively . 



Clearly, implementing the EFT approach to the AS = HWI requires carrying 
out new experiments in few-body systems for which ab initio structure computations can 
be performed. As outlined in Reference (0), a program of such measurements exists in 
principle. From a practical standpoint, carrying it out will involve meeting a number of 
experimental challenges. In light of these new theoretical developments and experimental 
opportunities, we believe it is time to review the field of hadronic PV anew. Comprehen- 
sive reviews of the subject have appeared over the years, including the influential Annual 
Reviews article by Aldelberger and Haxton completed two decades ago (0). In what 
follows, we hope to provide the "next generation" successor to that work, updating the 
authors' analysis in light of new theoretical and experimental progress. Since our focus 
will be on new developments, we touch only lightly on older work that has been reviewed 
in Reference (jij) and elsewhere (fiof). Before doing so in detail, however, we find it useful 
to summarize the primary developments and shifts in emphasis that have occurred since 
Reference (jjj appeared: 

• The extensive development of xPT and NN EFT, together with substantial progress 
in performing lattice QCD simulations, has revolutionized our approach to treating 
hadronic physics. While the use of hadronic models can provide important physical 
insights, the present day "holy grail" is to derive first-principles QCD predictions 
for hadronic phenomena. At the time of the Adelberger and Haxton review, the 
quark model was still in vogue, whereas lattice QCD and hadronic EFTs had yet 
to realize their potential. Today, the situation is reversed. Indeed, in the case of 
AS = HWI, the use of a meson-exchange model for the NN interaction that 
entails a truncation of the QCD spectrum and contains effective couplings that 
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likely parameterize more physics than the elementary meson-nucleon PV interaction 
(e.g., 27r-exchange) obscures rather than clarifies the connection with the SM. We 
now know how to do better. 

• New experimental and technological developments have opened the way to perform- 
ing PV experiments in few-body systems. The landscape now differs substantially 
from that of the 1980 's, at which time it appeared that measuring a number of 
C(10~ 7 ) effects in few-body processes was impractical. Indeed, two decades ago, 
the presence of the nuclear enhancement factors made experiments with many- 
body nuclei such as 18 F more attractive than those in few-body systems. Since 
then, precise new measurements of 10~ 7 PV observables in pp scattering, not spin 
rotation, and polarized neutron capture on hydrogen have either been completed or 
are in progress, and plans are being developed for other similarly precise few-body 
measurements at NIST, LANSCE, the SNS, and IASA (Athens). As we discuss 
below, completion of a comprehensive program of few-body measurements is now 
a realistic prospect. 

• Enormous progress has been made in performing precise, ab initio calculations in 
the few-body system using Green's function and variational Monte Carlo methods. 
These computations start with state-of-the-art phenomenological potentials that 
incorporate our present knowledge of NN phase shifts and include minimal three- 
body forces as needed to reproduce the triton binding energy and other three-body 
effects. A marriage between the NN EFT methods and these few-body computa- 
tional approaches is also being developed. As a result, a realistic prospect exists 
for performing precise computations with the PV EFT for few-body observables, 
leaving one free from the nuclear structure questions that enter the interpretation 
of many-body PV observables. 

In short, the frontier today for understanding the AS* = HWI lies in the few-body 
arena, for which a combination of precise experiments and first-principles theory provide 
new tools for making the most direct possible confrontation with the interplay of the 
strong and electroweak sectors of the SM. In the remainder of this article, we elaborate 
on this view. 



2 Weak Meson Exchange Model Meets the End of 
the Road 

While the era of the meson-exchange framework for hadronic PV is drawing to a close, 
it has played such a central role in the field that its development and use following the 
publication of Reference (0) calls for a brief review. The primary theoretical develop- 
ments have included updated theoretical "reasonable ranges" and "best values" for the 
h\ 4 provided by DDH and others (0, [ll| the analysis of nuclear anapole moments 
extracted from atomic PV experiments, computations of nuclear PV contributions to PV 
electron scattering asymmetries, and new global fits of the h l M to nuclear and hadronic 
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PV data. Experimentally, one has seen the completion of the TRIUMF 221 MeV pp 
scattering experiment and a neutron spin rotation experiment at NIST, the launching of 
an np — > d'y experiment at LANCSE, and the first non-zero result for a nuclear anapole 
moment in an atomic PV experiment with 133 Cs. 



2.1 Meson Exchange Model of the Weak N-N Interaction 

The meson-exchange, PV NN potential, V^h, is generated by the meson-exchange dia- 
grams of Figured, wherein one meson-nucleon vertex is parity conserving and the other 
parity violating. The Lagrangians for each set of interactions have been written down on 
numerous occasions in the literature, so we only give the final form of the static potential: 
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Here Pi = — iVi, with Vj denoting the gradient with respect to the coordinate Xi of the 
i-th nucleon, r = \x\ — x%\ is the separation between the two nucleons, 
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is the standard Yukawa function, and the strong ttNN coupling g n NN has been expressed 
in terms of the axial-current coupling qa using the Goldberger-Treiman relation: g n NN = 
gAiriN/F n , with F n = 92.4 MeV being the pion decay constant. The gy, V = p, u>, are 
the strong vector meson-nucleon Dirac couplings, and the xv give the ratio of the strong 
Pauli and Dirac couplings. The terms in Eq. (J2J) display different dependences on isospin 
and spin, so that various observables are sensitive to distinct linear combinations of the 
h l M . A notable feature is the absence of a neutral 7r-exchange component. Indeed, the 
only manifestation of 7r-exchange appears in the first term of Eq. (0) that contains only 
products of the isospin raising and lowering operators for the two nucleons. This feature 
reflects a more general theorem by Barton that forbids a neutral pseudoscalar-exchange 
component in the PV potential when CP is conserved ffl3h . 
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Table 1: Theoretical reasonable ranges (second column) and best values (columns 3-5) for 
the PV meson-nucleon couplings fla) . h % M , from DDH Dubovic and Zenkin (DZ) \tk) . 
and Feldman et al. All values are quoted in units of = 3.8 x 10~ 8 . 



PV Coupling DDH Range DDH Best Value DZ FCDH 



K 


- 


> 30 


+ 12 


+3 


+7 


h° 


30 - 


* -81 


-30 


-22 


-10 


K 


-1 - 


-> 


-0.5 


+ 1 


-1 


K 


-20 - 


-> -29 


-25 


-18 


-18 


hi 


15 - 


* -27 


-5 


-10 


-13 


hi 


-5 - 


-> -2 


-3 


-6 


-6 



The values of the h l M appearing in are most conveniently expressed in units 

of g n , the natural strength for the weak AS = 1 B — > B'tt couplings 2 : 

^ = 3.8xlO- 8 ^^S . (4) 

The original DDH reasonable ranges and updated best values are given in Table ^ Note 
that no prediction for hi' appears, as DDH were unable to compute this constant in 



Reference (p). Subsequently, Holstein (|14h used a 1/2 pole model to estimate this 



parameter. Using the quark model to compute the l/2~ <-> l/2 + mixing matrix elements, 
he obtained h 1 ' ~ 1.8 g n . Henceforth, we will not refer to this prediction when referring 
to the DDH values. 

The various SU^)™ symmetry arguments, current algebra techniques, and quark 
model estimates that lead to the values in Tabled have been discussed in detail elsewhere 
(0 0), and since our emphasis lies on a new formulation in which these couplings do 
not appear, we do not revisit those discussions here. Instead, we concentrate on new 
applications of this framework. 

Anapole Effects 

Two particularly novel uses of the PV meson-exchange framework have been in the 
analysis of atomic PV experiments and PV electron scattering. Shortly after PV was 
observed in /z-decay and /3-decay, Zeldovich and Vaks pointed out that weak interactions 
could also induce a PV coupling of the photon and fermion (jlfih . Electromagnetic (EM) 
gauge invariance implies that the lowest dimension effective operator for this coupling 
has the form |l7l ) 

4 J v 7 = |^75^/^^ , (5) 
where F^ v is the EM field strength tensor, Fa is the anapole coupling, and A is an 



2 Here, B and B' denote octet baryons. 
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appropriate mass scale. This effective operator leads to the momentum-space interaction 



M e JJ = -^f u{p) [q 2 lfl - 4q,\ T5«(p)e"(?) , (6) 

where q = p' — p is the momentum of a photon with polarization vector e^{q). From (JUjK 
it is clear that the anapole coupling involves only virtual photons, since for a real photon 
q 2 = and one can always choose a gauge in which q ■ e = 0. Since d u F^ u = Eq. (0) 
implies that the anapole interaction effectively couples the fermion axial current to the 
source of the EM field, J M . 

In the 1980's, Flambaum, Khriplovich, and Sushkov (flii Hi! ) observed that the 
anapole moments of nuclei would scale as the square of the nuclear radius, rather than 
as 1/A 2 , so that their magnitudes would be enhanced as A 2 ^ 3 in heavy nuclei. Moreover, 
the nuclear anapole moment would couple the nuclear axial current to the EM currents 
of the atomic electrons, thereby inducing a PV, nuclear spin-dependent (NSD) term in 
the atomic Hamiltonian. Experimentally, one could isolate this effect by observing NSD 
transitions in atomic PV processes. As discussed below, a non-zero result for the 133 Cs 
anapole moment has been obtained by the Boulder group, while limits have been placed 
on the anapole moment of 205 T1 by the Seattle group. Efforts are presently underway to 
measure the anapole moments of other nuclei, such as francium (for recent reviews, see 
e.g. references Q 0)). 



F A 



The new anapole moment measurements have stimulated considerable theoretical 
activity. Using a one-body averaged version of V^ n and a simple single particle shell 
model, the authors of references llih estimated the magnitude of the anapole mo- 
ments of various nuclei, demonstrating the A 2 !' 3 scaling under these conditions. Substan- 
tially more sophisticated shell model calculations - using the complete two-body potential 
and associated meson-exchange currents - were carried out in references (|7|,|23j,|24). The 



results have been used to extract constraints on the DDH couplings, as shown in Figure 
12] • Of particular interest is the significant disagreement between the 133 Cs anapole con- 
straints on the relevant isovector combination of couplings, compared with those obtained 
from the circular polarization P 7 in the 7-decay of 18 F. Other nuclear model computa- 
tions of the 133 Cs anapole moment, havin g gr eater or lesser degrees of sophistication, lead 
to similar conclusions (@, H 0, HS B . 



As discussed in references (p, l24f ). state-of-the art shell model computations un- 
avoidably entail model-space truncations, and in the case of the cesium anapole moment 
calculations, inclusion of the omitted contributions would likely increase, rather than 
decrease, the disagreement with the 18 F result. In contrast, the nuclear structure analy- 
sis used to interpret the P 7 results is thought to be robust, since the dominant nuclear 



mixing matrix element can be calibrated against an analog /3-decay amplitude (|31l . 13^ 
Thus, the implications of the new result for the Cs anapole moment are quite puzzling. 

The AS = HWI can also contribute to the PV asymmetry Apy measured in the 
scattering of longitudinally polarized electrons from hadronic and nuclear targets. In the 
late 1990's, it was realized that the anapole moment of the proton contributes to Apy 
for elastic ep scattering in a way that is indistinguishable from that of the axial vector 
coupling of the Z° to the proton, or G\ (see Figure EJ) Moreover, it was shown - 

using the DDH framework - that the both the magnitude of the proton anapole moment 
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Figure 2: Constraints on effective DDH weak meson-nucleon couplings deduced from PV 
observables in nuclei and anapole moments of heavy atoms (courtesy of W. C. Haxton). 
Here, h\f c = h\ - 0.12/iJ - 0.18/^ and hf c = -(h° p + 0.7h°J 



contribution as well as the theoretical uncertainty associated with it was sufficiently 
large as to significantly affect the interpretation of A-py. At that time, a program of PV 
electron scattering measurements was being developed to determine the strange quark 
contributions to the electric and magnetic form factors of the proton, G S E and G S M , 
respectively. The presence of the anapole related uncertainties would be particularly 
problematic for the extraction of G S M from the backward angle asymmetry measurements 



Consequently, additional measurements of Apy for quasielastic (QE) scattering from 
the deuteron were carried out. Since the deuterium asymmetry is strongly sensitive to G e A 
but considerably less sensitive to G S M than is the proton asymmetry, a measurement of 
Apy(eD) - in conjunction with Ap\(ep) - could be used to test the theoretical estimates 
of Reference (js^h while providing for a determination of G S M that is independent of 
hadronic PV uncertainties. The initial results of these measurements, completed by the 
SAMPLE Collaboration (jH), yielded a new puzzle: the effective G e A extracted from the 
ep and eD asymmetries was consistent with zero. The calculations of Reference (Jiih had 
predicted a ~ 40% reduction to the value of G e A arising from SM electroweak radiative 
corrections and the anapole effect. The SAMPLE result, however, implied a substantial 
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enhancement of the anapole contribution or that of other radiative corrections over the 
predictions of Reference (jiih . 




Figure 3: Contributions from hadronic PV (proton anapole moment) to the effective axial 
vector electron-proton coupling, G e A . 



Subsequent theoretical studies attempted to determine the origin of the anomaly, 
scrutinizing various contributions to Ap\: the original computation of Reference (Bs ) was 
revisited and updated using heavy baryon %PT (1371 ): possible quark model enhancements 
were considered (|38f ): the g 2 -dependence of the anapole contribution was studied f 39l 1401) . 
and contributions from parity mixing in the deuteron and final np states were computed 
using the meson-exchange model (j4ll |42J). In all cases, no large effects were found that 
could resolve the puzzle. Ultimately, the SAMPLE Collaboration carried out a reanalysis 
of the pion-production background in the deuterium experiment that shifted the value 
and brought the axial term into agreement with the predictions of 
The theoretical results from the latter work have now been used in 



of A^(eD) 
references 



>Lj43|) a 
(33, 37). 



extracting G S M from the backward angle proton asymmetry (j44l ). 

Within the DDH framework, these developments - along with the completion of a 
new pp scattering experiment at TRIUMF - have motivated reanalyses of the hadronic 
and nuclear PV observables in terms of the h l M . After discussing the recent experimental 
developments, we provide a summary of our current understanding of hadronic PV in 
this context and make the case that a fundamental paradigm shift is required in order to 
make further progress in this field. 

2.2 Experimental Progress 

Earlier reviews, e.g. (0), |l5h. (|45h. have documented an extensive body of experimental 
work aimed at characterizing the AS = HWI, largely carried out in many-body nuclei. 
As noted in the Introduction, the PV effects in the much simpler NN and few nucleon 
systems were almost unpractically small [(9(10 -7 )] from the standpoint of past experi- 
mental feasibility, so the realization of fortuitous nuclear structure effects that 'amplify' 
the underlying NN PV signal by several orders of magnitude naturally led to an earlier 
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focus on many-body systems. Future progress will hinge on a handful of precise exper- 
iments in much simpler few nucleon systems, for which the theoretical interpretation 
is less fraught with mo del- dependent uncertainties and for which precise measurements 
now appear to be realistic. Here, we review recent experimental progress in the few-body 
sector and comment on new developments in probes of PV in many-body nuclei. These 
developments have been cast largely in the weak meson-exchange framework, so we will 
use it here in discussion their theoretical implications. When developing the EFT frame- 
work in Section |31 we will reframe the discussion of these theoretical implications in the 
EFT formulation. 

Probing the NN weak interaction in few body systems presents significant exper- 
imental challenges. The bare NN O(10~ 7 ) PV effects compete with a host of potential 
systematic errors at this level, which must be both minimized through careful experimen- 
tal design and simultaneously measured to ensure that they do not obscure the true PV 
signal. Even acquiring sufficient data to reach a statistical error at the 10~' — 1CT 8 level 
is no mean feat, necessitating the use of current mode detection which in turn introduces 
its own systematic error sensitivities that must be controlled and understood. The most 
accessible NN and few nucleon observables are accessed in polarized beam experiments, 
where rapid polarization reversal provides a practical means of suppressing low frequency 
noise and systematic effects, at the expense of introducing a built in sensitivity to spin 
correlated beam properties that can mimic the PV observable in question. In such cases, 
as much effort must be expended to optimize and characterize the polarized beam prop- 
erties as is required to design and commission the PV experimental apparatus. Typically, 
a successful experiment spans a decade or more from initial concept to publication of a 
significant result, with continuous refinements of the experimental apparatus and tech- 
nique until the desired sensitivity is reached. In most cases, an order of magnitude more 
data are required to refine and test the apparatus than to acquire the final PV data 
sample. 

To date, there have been a number of significant measurements of PV in pp scat- 
tering, but despite several decades of experimental effort, a definitive observation of PV 
in the np system remains to be established. Recent experimental progress in the NN 
system includes the completion of a program of high precision measurements of PV in pp 
scattering which yield independent constraints on the weak couplings of heavier mesons, 
and commissioning of a PV asymmetry measurement in np — > d^j that is aimed at a 
precise determination of h\. These experiments will be discussed briefly below; future 
possibilities for precise measurements in two and few-nucleon systems are discussed in 
Section 3. 

2.2.1 Longitudinal Analyzing Power in pp Scattering 

The PV observable that has been studied in pp scattering is the longitudinal analyz- 
ing power, A z = , where a + and o~ are the elastic scattering cross sections for 
positive and negative helicity beams incident on an unpolarized hydrogen target. The 
analyzing power can be naturally expressed as a sum of parity mixed partial wave contri- 
butions, with only S — P mixing required to characterize A z at low energy. The first two 
partial waves in the expansion are sufficient to describe A z up to several hundred MeV; 
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the ( 1 5o — 3 Pq) contribution dominates at low energy, while the ( 3 P2 — 1 Dq) amplitude 
starts to become significant above 100 MeV. These two partial wave amplitudes have 
complementary dependences on the weak M-N couplings: h p p p = [h° p + hi + h 2 p /\/6) and 
h v J = h%j + /i*. As noted by Simonius (46), the dominant S — P wave mixing integrates 
to zero near 220 MeV beam energy, due to a fortuitous cancellation of the strong S 
and P wave phase shifts. This observation motivated the recently completed TRIUMF 
experiment at 221 MeV (|47l). which was designed to isolate ( 3 P2 — 1 D2) contribution. 
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Figure 4: Principle of the TRIUMF p-p parity violation experiment. Longitudinally po- 
larized protons at 221 MeV passed through a 40 cm liquid hydrogen target, which scattered 
~ 4% of the beam. The polarization-dependent target transmission was measured by per- 
forming an analog subtraction of two dc current signals from transverse field ionization 
chambers. Reprinted with permission from Berdoz et at, Phys. Rev. C Vol. 68, 034004 
(2003), Fig. 1. Copyright (2003) by the American Physical Society. 



While earlier high precision measurements at low energy were performed in a total 
scattering geometry, the TRIUMF 221 MeV measurements were carried out in trans- 
mission mode. The small size of the total scattering asymmetry A z ~ 10~ 7 implied a 
transmission asymmetry of order 10 -9 for the TRIUMF experiment. Integrating detec- 
tors with small angular acceptance coupled to low noise electronics, excellent beam and 
liquid hydrogen target stability, and a highly polarized beam with minimal helicity cor- 
related beam properties were essential to the success of the measurements. Rapid beam 
polarization reversal (40 Hz, with controlled phase slip with respect to the line frequency) 
led to an ac parity violating signal at a well determined frequency, greatly suppressing 
the noise contributions with respect to a dc measurement (which would be impossible in 
this case due to the extremely small size of the parity violating asymmetry signal). The 
principle of the TRIUMF transmission mode measurement is illustrated in Figure EJ 
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The TRIUMF laboratory spent many years developing a state-of-the-art optically 
pumped polarized H- ion source (OPPIS) (|4a. Eifl. which was ideal for demanding sym- 
metry tests that require polarized beam, since the polarization (~ 85%) was reversed 
by changing the frequency of the laser light with no changes to macroscopic electric or 
magnetic fields that could influence the beam properties. Even so, significant correc- 
tions (~ 40 %) had to be made to the raw asymmetry for helicity correlated transverse 
polarization components. 

Transverse polarization is a pathological source of systematic error in pp scattering 
measurements because of the relatively enormous parity-allowed transverse analyzing 
power A y , which is over a million times larger than the PV longitudinal analyzing power 
A z . While it is not practical to demonstrate and maintain tiny transverse polarization 
components at or below the 10~ 6 level, fortunately the geometrical symmetry of the 
apparatus can be invoked to establish a 'neutral axis' for the beam transport such that 
the false asymmetry arising from transverse polarization components is identically zero 
if the beam is locked on this axis by a position feedback system. This technique was 
used to ensure that corrections for average transverse polarization components (P y ) and 
(P x ) were insignificant at the 10~ 9 level in the TRIUMF measurements - the associated 
corrections AA Z e.g. for transverse vertical polarization (P y ) scale as (x)(P y ), where 
(x) is the net displacement from the symmetry axis, whose location can be determined 
by separate calibration experiments. Note that the form of the driving term (x)(P y ) 
is that of the first moment of transverse polarization at the detector location, which 
is referred to as an extrinsic first moment; a similar correction must also be made for 
intrinsic first moments of the form (x P y ) which arise from nonuniform distributions of 
transverse polarization within the beam envelope. These extrinsic polarization moments 
are extremely difficult to measure and control. 
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Figure 5: Demonstration of a new current mode scanning polarimeter at TRIUMF. Left 
panel: intrinsic polarization moments (xP y ) measured simultaneously with counting mode 
(vertical axis) and current mode (horizontal axis) devices. The errors on the current mode 
measurements are too small to display on this scale. The dotted line has a slope of unity, 
and the solid line is a fit to the data. Right panel: instrinsic polarization moments for a 
given beamline tune, measured with a pair of current mode PPMs. The ratio of moments 
at PPM1 and PPM2 remains relatively constant, thus verifying an essential assumption 
of the systematic error reduction scheme. Reprinted with permission from Ramsay et al., 
Proc. 9th International Workshop on Polarized Sources and Targets, p. 289-293 (2002), 
Fig. 3. Copyright (2002) by World Scientific. 
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In the TRIUMF experiment, a pair of scanning polarimeters (|50l ) was employed 
to provide continuous measurements of the distribution of transverse polarization com- 
ponents within the beam, interleaved with parity data taking from the transverse ion 
chambers on an 8-state, 200 ms data cycle. The limiting factor in the TRIUMF ex- 
periment was the statistical precision of the PPMs, which were coincidence counting 
mode devices. An essential assumption made in determining the correction for intrinsic 
polarization moments was that the moments evolved linearly with position along the 
beamline and thus had a stable ratio between upstream and downstream PPMs. This 
ratio could in principle be tuned to achieve a null sensitivity to intrinsic polarization mo- 
ments, an approach modelled on the successful polarization neutral axis idea described 
above. Unfortunately, the ability of the PPMs to measure this first moment ratio suffi- 
ciently precisely in a reasonable amount of time was severely limited, and so significant 
corrections still had to be made to the data. 

R&D efforts towards a follow up experiment at higher energy (which was never 
realized) included development of a current mode scanning polarimeter which had ~ 
20 x greater statistical precision (pll); the current mode polarimeter was able to clearly 
demonstrate the constant linear evolution of intrinsic moments along the beamline, as 
illustrated in figure thus independently validating the corrections procedure that was 
used to obtain a final result for A z at 221 MeV from TRIUMF experiment 497. 
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Figure 6: The most precise measurements of parity violation in pp scattering at low and 
intermediate energy, and recent theoretical predictions. Experiments we re p erformed at 
Bonn (13.6 MeV) fM), PSI (45 MeV) @j and TRIUMF (221 MeV) (Qj. The solid 
curve shows the calculation by Carlson et al. 1^) including a fit of the weak meson- 
nucleon coupling constants to the data. Reprinted with permission from Berdoz et al., 
Phys. Rev. C Vol. 68, 034004 (2003), Fig. 13. Copyright (2003) by the American 
Physical Society. 
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scattering has attracted considerable theoretical interest since 
. Recent calculations are shown in Figure El together with the 
recise experimental data at low and intermediate energies. Driscoll and Miller 
used the Bonn potential to treat the strong NN interaction, with weak meson- 
nucleon couplings taken from Reference (0). Iqbal and Niskanen's calculation adds a A 
isobar contribution (57) to the Driscoll and Miller model. The calculation of Driscoll and 
Meissner (Joli ) is based on a chiral soliton model, while the quark model calculation of 
Grach andThmatikov © takes explicit account of quark de grees of freedom. 

Figure [3 by Carlson et al. p4j ) shows the limits on the weak meson-nucleon cou- 
plings h pp and h vp imposed by the low energy pp asymmetry measurements and the 221 
MeV TRIUMF result. The error bands are based on a calculation assuming the Argonne 
fis (AV-18) potential (j^), the Bonn 2000 (CD-Bonn) (joll ) strong interaction coupling 
constants, and including all partial waves up to J=8. Although the TRIUMF measure- 
ment is not sensitive to A z from SP mixing, and the contribution from PD mixing 
contains no h%f contribution, there is some h vp dependence arising from higher partial 
waves. The best fit to the pp data yields hf = -22.3 x 10~ 7 and h™ = 5.17 x 10~ 7 , 
compared to the DDH "best guess" values of h pp = -15.5 x 10" 7 and hff = -3.0 x 10~ 7 . 




Figure 7: Present constraints on the weak meson-nucleon couplings h pp and h pp , fitted to 
the low energy and 221 MeV pp asymmetry data ( fj^) . The plot shows curves of constant 
total x 2 = 1,^,3,4 an d 5. Axis scales are 10~ 7 . Reprinted with permission from Carlson 
et al, Phys. Rev. C Vol. 65, 035502 (2002), Fig. 8. Copyright (2003) by the American 
Physical Society. 
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2.2.2 Progress in the np System 

Parity violation in the np system can in principle be detected in a variety of processes 
that can reveal complementary aspects of the weak NN interaction. In np capture, there 
are two complementary PV observables: P^, the 2.2 MeV 7 ray circular polarization for 
an unpolarized neutron beam, and At, the asymmetry in the emission of 7 rays with 
respect to the neutron spin direction if the beam is polarized. Closely related to the first 
of these, P^, is the helicity asymmetry A\ in the photodisintegration of deuterium with 
circularly polarized photons; the two are asymptotically equal to each other at threshold, 
while A\ is predicted to drop rapidly with increasing photon energy, fal ling an order 
of magnitude as the photon energy increases to 1 MeV above threshold (|62l). Finally, 
the transmission of polarized neutrons through hydrogen should reveal a tiny PV spin 
rotation about the neutron propagation direction z: d<f) np / 'dz . 

Unfortunately, all of these ny s yste m measurements are extremely challenging; the 
first three have been attempted (0, but have yielded null results with limits at 
least one order of magnitude too lar ge to provide a meaningful constraint on the weak 
meson-exchange predictions (|62l . I66L1671 loa . l69f ) all of which are at the level 5 x 1CT 8 
or smaller. The two PV observables involving the np capture reaction, P 7 and have 
complementary dependences on the weak MN couplings - notably, can yield a unique 
constraint on h\, while P^ depends on a linear combination of ir, p and u weak couplings 
(0). Measurement of gamma ray circular polarization requires a Compton polarimeter 
with typical sensitivity at the few percent level, thus rendering the P^ measurement 
in principle even more challenging than a measurement of A^. A measurement of the 
helicity asymmetry A\ in the photodisintegration of deuterium is at the early conceptual 
design stage for the future experimental program at IASA, Athens fz3), with the aim 
of reaching a sensitivity at the 10~ 8 level. The np spin rotation measurement, with an 
anticipated d<p np /dz = 5 x 10~ 7 rad/m (1711 ) has not yet been attempted, but is envisioned 
as a future component of the SNS fundamental neutron physics program, as discussed in 
Section 13.41 below. 

The intrinsic interest of a clean measurement of h\ in the NN system, together with 
considerations of experimental feasibility, have led to the launching of a major effort at 
LANSCE 0) - the NPDGamma experiment - to make a definitive measurement of A , 
with an ultimate goal of reaching ± 10% of the DDH prediction. The measurement of A^ 
requires a polarized neutron beam with precisely known spin direction and a measurement 
of the 2.22 MeV gamma ray angular distribution: % ~ (1 + A~ t cosO), where 9 is the 
angle between the gamma ray momentum and the neutron spin. Even at milli-eV neutron 
energies for neutrons moderated in liquid hydrogen, the scattering cross section exceeds 
the capture cross section by a significant factor. A parahydrogen target is essential to 
avoid neutron spin flip on scattering in the target, and the useful beam flux is below 15 
meV to prevent depolarization in the target. 

A previous measurement of A^, performed at the ILL reactor in the 1970's 
reported a value of (0.6 ±2.1) x 10~ 7 ; the NPDGamma experiment is designed to reach 
an ultimate sensitivity of ±5 x 10 -9 , with the uncertainty dominated by statistical rather 
than systematic errors. Major improvements in the experimental instrumentation and 
techniques that should make this possible include: 



18 



• use of a high intensity, low energy, pulsed beam, which allows for neutron energy 
determination via time-of-flight measurement and allows for a separation in time 
of prompt 7 ray background from the neutrons of interest; 

• polarization of the beam via selective transmission through optically pumped po- 
larized cells of 3 He; the well known spin- dependent cross section leads to an energy 
dependent polarization that can be very high for low energy neutrons, and the beam 
polarization can be continuously and directly monitored by online measurements 
of the 3 He cell transmission; 

• use of a resonant RF spin flipper capable of flipping spins at all neutron energies 
with high efficiency and eliminating Stern-Gerlach steering of the neutron beam 
associated with spin flip; 

• implementation of a large solid angle, high efficiency Csl (Tl) gamma detector array 
(|73i ) instrumented with sensitive current mode electronics whose intrinsic noise is 
negligible compared to neutron counting statistics. 
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Figure 8: Layout of the NPDGamma apparatus currently taking data on Flight Path 12 
at LANSCE. A highly uniform, vertical magnetic guide field to preserve the neutron spin 
direction is provided by a set of field coils (not shown). 

A schematic of the NPDGamma apparatus mounted in the new experimental cave 
on flight path 12 at LANSCE is shown in Figure |H1 The neutron spins, polarized by the 
3 He transmission cell, are efficiently reversed by the RF spin flipper on a pulse-by-pulse 
basis, thus alternating the sign of the parity violating asymmetry measured in the gamma 
detector array at 20 Hz. At the time of writing, NPDGamma has achieved a number of 
major milestones and is ready to take production data at LANSCE, pending installation 
of the liquid hydrogen target. The apparatus has been fully commissioned; electronic 
asymmetries are consistent with zero at the few x 10 -9 level, and PV asymmmetries 
arising from neutron capture on a variety of solid targets which are representative of 
materials used to construct the beamline and parity instrumentation have been measured. 
All are consistent with zero at the 10~ 6 level or smaller, with sufficient accuracy to 
conclude that background asymmetries will not play a significant role in the hydrogen 
target data. A known PV asymmetry in CI at the 10~ 5 level has been remeasured and 



19 



will be used as a diagnostic tool for NPD Gamma, with the periodic insertion of a CCI4 
target to verify the consistent performance of the experimental setup. Contamination of 
the hydrogen target PV (up-down) asymmetry by the comparably small, parity-allowed 
left-right asymmetry will be kept below an acceptable level by determining the effective 
detector angles in situ - this in turn will be accomplished by scanning the detector array 
horizontally and vertically with respect to the target while acquiring np capture data. 
The measured chlorine asymmetry is shown in figure the net PV (up-down) and parity- 
allowed left-right asymmetries deduced from these data ((73) are: A 1 = (—19 ± 2) x 1CT 6 
and A LR = (-1 ± 2) x 10" 6 . 



x1Q 3 




Detector Angles (Averaged Over 4 Rings) [radians] 

Figure 9: CCI4 gamma ray asymmetries, calculated from opposing detector pairs, plot- 
ted versus angle of the first detector in the pair with respect to the vertical, for the 
NPDGamma 2004 data. The total asymmetry A = A^cosO + AiRsinO is deduced from 
the fit. Figure courtesy of M. T. Gericke. 

Unfortunately, a number of factors have conspired that severely limit the statistical 
accuracy that can be achieved by running the NPDGamma experiment at LANSCE. 
These include a factor of 4 reduction in the available neutron flux per beam pulse as 
compared with expectations prior to the recent upgrade of the LANSCE facility, and a 
further reduction in high quality beam time due to magnetic interference from a neigh- 
boring experiment. In view of these limitations, the NPDGamma collaboration plans 
to carry out a first measurement with hydrogen in 2005-2006, which would provide a 
statistics limited result for A d accurate to ±1 x 10~ 7 or better. At the time of writing, 
plans are being made to move the experiment to the Fundamental Neutron Physics Beam 
Line at the Spallation Neutron Source, which would enable the collaboration to make a 
measurement of A d to ±1 x 10 -8 or better, as discussed in Section EH 
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2.2.3 Neutron Spin Rotation Experiments 

Low energy neutrons can exhibit a parity violating spin rotation induced by the hadronic 
weak interaction as they pass through a medium. The observable, 4>pv, is the angle 
of transverse spin rotation about the neutron's direction of motion. In the limit of 
zero neutron energy, the PV neutron spin rotation is energy independent (jzij). After 
traversing a distance z through the medium, the neutron spin will precess by an amount 
4>pv = Znpzfpy, where fpy is the parity violating coherent forward scattering amplitude 
for a low energy neutron and p is the number density of the medium ([Til ). The basic 
requirements of a neutron spin rotation experiment are a source of low energy polarized 
neutrons, a target containing the material of interest, and a spin analyzer downstream 
of the target, from which the value of (fipy can be deduced. 

Two cases are particularly interesting from the standpoint of testing models of the 
hadronic weak interaction: parity violating neutron spin rotations in ^-H and 4 He. These 
two cases have complementary dependences on the weak meson-nucleon couplings. The 
np case is dominated by the weak pion exchange contribution, and thus would yield 
similar information to the NPDGamma experiment described above. The 4 He PV spin 
rotation has a significant contribution from h\ but is also sensitive to the h° p coupling. 
Experimentally, the 4 He case is more tractable, due to a much longer mean free path for 
low energy neutrons than hydrogen. Like NPDGamma, the np spin rotation experiment 
would require a parahydrogen target to avoid neutron depolarization due to spin-exchange 
collisions with the hydrogen molecules, and only neutron energies below 15 meV would 
be useful for the experiment. 

A 4 He spin rotation measurement was carried out at NIST in the 1990's |7fih . with 
the result: <fipv = (8 ± I4: sta t ± 2 sys ) x 1CT 7 rad/m. Unfortunately this heroic effort 
did not reach sufficient sensivitity to test hadronic weak interaction models. A possible 
experiment in hydrogen has been considered but not yet proposed, and will hopefully be 
carried out at the SNS within the next decade (see Section 3.3). 

What confounds the beautiful simplicity of the measurement principle is the ex- 
tremely small size of the PV rotations expected based on the DDH model: for 4 He, the 



predicted effect is <ppv = — 1 x 10 7 rad/m (|77l). while for the proton, the effect should 



be about 5 times larger (|7l[). These PV spin rotations are 8 orders of magnitude smaller 
than those induced by the earth's magnetic field in a typical measurement apparatus! 
Therefore, a much more elaborate scheme is required in order to achieve the necessary 
sensitivity. Great care must be taken to reduce ambient magnetic fields by many orders 
of magnitude; residual fields must be monitored carefully, and the experiment should be 
designed so that the effects of residual magnetic fields are cancelled to the maximum 
possible extent in the extraction of the PV spin rotation angle. 

The basic experimental technique employed at NIST and foreseen for subsequent 
measurements involves a double beam / double cell apparatus, as sketched in Figure ITU1 
The double cell design incorporates a 180° spin precession about the initial spin direction 
between the two target cells, one of which is full and one empty, for a given measurement. 
A clear advantage of the double cell design is that the much larger spin rotations due 
to residual magnetic fields exactly cancel in the two target states, to the extent that 
the magnetic fields and the neutron trajectories are exactly the same in both cases. An 
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Figure 10: Schematic of the double beam / double target system for 4 He neutron spin 
rotation measurements. The incident neutron spins are polarized out of the page. The 
PV effect is a spin precession around the direction of motion, which leads to components 
in the horizontal plane as shown. Either the front or the back target is full (shaded), 
and the other is empty, in channels A and B. Much larger spin rotations due to local 
magnetic fields are independent of the target state and are cancelled by subtraction. The 
it coil situated between the front and rear targets precesses the neutron spin by 18CP about 
the initial spin direction, thus reversing the sign of the PV spin rotation relative to that 
due to magnetic fields, indicated in the figure; the counting rate (B - A) — 2 4>pv- (Figure 
courtesy of Anna Micherdzinska) 



identical double target system sits beside the first one in the same cryostat, half the beam 
passes through each, and the two systems are run so that at any given time, the sign of the 
PV spin rotation in the two subsystems is exactly opposite. Downstream of the double 
target system is a polarization analyzer and segmented 3 He ionization chamber which 
measures the counting rate in 'left' and 'right' elements of the apparatus simultaneously. 

Instead of rapidly flipping the beam polarization as in the pp and NPDGamma 
experiments described earlier, the sign of the PV effect is reversed by alternating the 
location of the full target cell with respect to the 7r coil. The emptying and filling of the 
target cells took several minutes in the first spin rotation experiment, and each measuring 
interval was 10 min, corresponding to a reversal rate of 1 mHz and roughly 25% dead 
time. However, with a pair of double cell setups located side by side in the same cryostat, 
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and a beam splitter upstream of the apparatus, reactor beam intensity fluctuations can 
be effectively cancelled, which makes up for the slow PV reversal frequency. 

Features of the upgraded apparatus include: improved magnetic shielding, resulting 
in a fourfold reduction of ambient fields to the 20 /zG level; improved cryogenic systems 
with a faster target cycling time; operation with superfluid 4 He, which has a much smaller 
total and small-angle neutron scattering cross section than normal 4 He, and improvements 
to the NIST cold source leading to a 50% increase in the beam fluence. These and other 
improvements should result in a statistical error of ±3 x 10~ 7 rad/m in a 3 month run at 
NIST during 2005-06, with systematic errors (based on simulations of the experiment) 
at the 10~ 8 level. Like NPDGamma, the collaboration proposes to move the 4 He spin 
rotation experiment to the SNS, where the high intensity pulsed cold neutron beam 
would permit a high statistics measurement in about a year's running time with greater 
diagnostic capability for systematic errors (fTcl). 




2.2.4 Nuclear Anapole Moments 

Nuclear anapole moments, which may be accessed via measurements of parity violation in 
atoms, have recently opened a new window on the hadronic weak interaction This 
subject has been the focus of significant theoretical and experimental effort, including a 
recent review article in this journal (0). A detailed discussion of experimental atomic 
physics approaches is beyond the scope of the present work; however, for completeness, 
we include a brief overview of experimental efforts here. 

A prime motivation for atomic PV measurements is to provide stringent tests of the 
electroweak Standard Model. The dominant contribution to the neutral weak electron - 
nucleus interaction arises from an axial coupling to the electron and a vector coupling to 
the nucleus and is proportional to the nuclear weak charge, Qw The value of Qw can 
be computed with high accuracy in the SM, and a precise experimental determination 
of the nuclear weak charge can probe for deviations that one might expect to arise from 
physics beyond the SM. This dominant PV contribution is nuclear spin independent, 
and the extraction of Qw is minimally affected by hadronic weak interaction corrections. 
However, the weak electron - nucleus interaction also has a contribution from the vector 
coupling of the electron to the axial current of the nucleus; this leads to a nuclear spin 
dependent (NSD) parity violating interaction that can be detected as a small hyperfine 
dependence of atomic PV. The parity violating interaction in this case has contributions 
from Z° exchange, a hyperfine interaction correction, and finally the nuclear anapole 
moment (0,113). 

The weak interaction typically induces very tiny opposite-parity wavefunction ad- 
mixtures in atomic states, on the order of 10~ n . These in turn give rise to extremely 
small and otherwise parity-forbidden transition amplitudes. The best explored case is 
atomic PV in Cs, measured in the highly forbidden 6s- 7s transition in an atomic beam, 
where the tiny PV amplitude was arranged to interfere with a much larger Stark-induced 
transition amplitude in the presence of a static E field. 3 The authors (|79T ) report the use 

3 The advantage of this seemingly complicated experimental scheme is that the interference technique 
yields an observable which is linear in the tiny PV amplitude, rather than trying to measure the forbidden 
transition rate, which is quadratic. 
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of 31 different servo systems to precisely control electrical, optical and mechanical sys- 
tems during the measurements; they designed an experiment with 5 independent means 
of reversing the parity signal, and spent about 20 times more data taking effort in the 
investigation and elimination of systematic errors compared to actual PV data taking. 
From the standpoint of atomic structure calculations needed to interpret the measure- 
ments, the heavy alkali atom Cs is ideal for precision PV studies, due to its relatively 
simple electronic configuration of one valence (6s) electron outside a tightly bound Xe 
noble gas core. 

Since the nuclear anapole contribution is expected to scale as A 2//3 , and the atomic 
PV signature used to detect it scales as Zs, heavy nuclei are the preferred systems for 
experimental studies. To date, a definitive experimental result exists only for the anapole 
moment of 133 Cs (Fih . while an upper limit has been obtained for the anapole moment of 
Tl (30% 203 T1, 70 % 205 T1 isotopic ratio) (0) obtained via PV optical rotation measure- 
ments. In Cs, the nuclear spin independent atomic PV effect has been measured to 0.4%, 
while the anapole moment contribution has been determined to 14%. In Tl, the spin in- 
dependent effect has been measured to 3%, and the spin dependent anapole contribution 
is consistent with zero. These are the result of truly heroic experimental investigations 
spanning over a decade of instrumentation development and testing. Interpretation of 
the spin independent PV as a Standard Model test requires excellent understanding of 
the atomic structure, while the anapole moment must be first unravelled from the com- 
peting effects of Z° exchange and a hyperfine correction, and then interpreted in terms 
of a HWI model, e.g. the meson-exchange model of DDH. 

Atomic PV has been measured in other systems, but to lower precision than in 
Cs, and with greater uncertainty in the atomic structure details needed to interpret the 
measurements. Currently, a new generation of experiments is under development, aimed 
at determining both the nuclear spin independent and anapole moment contributions 
to high precision. This includes ongoing efforts to develop alternative experimental ap- 
proaches to atomic PV measurements in Cs; demonstration of a new technique based on 
stimulated emission has recently been reported (1811) . anticipating that an ultimate pre- 
cision of 0.1% could be reached. As remarked in (|82l). there is a significant advantage to 
measuring PV in a range of isotopes of the same atom, since a number of atomic structure 
uncertainties cancel when PV ratios are considered. Atomic ytterbium (Z=70), having 
7 stable isotopes, has been proposed (Jiih , and preparatory spectroscopic studies are 
underway (|84l). Besides the atomic beam techniques described earlier, new experimental 
approaches based on trapped atoms or ions are also being explored (85, 86). 

Perhaps the 'holy grail' of this field is the study of atomic PV in radioactive francium 
isotopes (Z=87), the heaviest alkali atomic system, recently reviewed in Reference (jl3). 
Both the spin independent and anapole PV contributions are expected to be roughly an 
order of magnitude larger for Fr (depending on the isotope) than for Cs. In addition, 
Fr has a large number of isotopes spanning almost 30 neutrons with lifetimes greater 
than 1 s that cover a wide range of nuclear structure conditions, which in principle 
permits an unprecedented systematic study of atomic PV in a simple atomic system. 
With the longest lived isotope 223 Fr having a half life of only 23 minutes, a radioactive 
beam facility is required for a future program of Fr PV studies. Pioneering atomic 
spectroscopy studies of a number of trapped Francium isotopes have been carried out at 
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Stony Brook (|88l). and a collaboration has recently formed with the goal of establishing 
a long term francium program at the TRIUMF ISAC radioactive beam facility (89), 
where an actinide production target is planned and necessary to produce the required 
quantities of Francium (and other heavy systems that are of interest for fundamental 
symmetry tests such as proposed EDM measurements in radon). 



2.2.5 Parity Violation in Compound Nuclei 

In contrast to the extremely small (~ 1CT 7 ) PV asymmetries in the two nucleon system, 
many-body nuclei have over the years provided many examples of parity violation at a 
much more significant scale, some of them surprisingly large (see e.g. Reference (Jjj) and 
references therein). For the handful of cases for which reliable nuclear wavefunctions can 
be used to interpret the data, there is qualitative agreement with predictions of the meson- 
exchange model as illustrated in Figure 121 albeit with a significant discrepancy in the 
scale of the pion coupling as discussed earlier. Since the publication of Reference (0), an 
extensive and in many ways complementary program of PV measurements in compound 
nuclei has been carried out, dominated by the work of the TRIPLE collaboration at 
LANSCE. A review of this field was presented in this journal in 1993 (J45J) focusing on 
resonances in 238 U and 232 Th; an updated c omp rehensive review of that and more recent 
data and analysis is provided in Reference (90). 

The TRIPLE collaboration measured longitudinal cross section asymmetries for 
neutron energies in the eV - keV range on a range of nuclear targets. Mixing of S- 
and P- wave resonances of the same J leads to large PV effects in some cases, with 
enhancement factors as large as 10 6 relative to the nucleon-nucleon PV asymmetries, due 
to the high density of states in the compound nucleus. The measurement of a number of 
PV asymmetries in the same nucleus is a key feature of this program, which is essential 
to the interpretation of the results. Using a statistical approach, the collaboration was 
able to extract results for either the RMS weak mixing matrix element Mj or the weak 
spreading width r<„ = 2iiMj/Dj where Dj is the level spacing, based on a total of 75 
PV resonances observed in 18 nuclei. Important spectroscopic information (J values) for 
the resonant states was known in some cases, and modelled in other cases, in order to 
interpret the results. The weak spreading widths were found to be on the order of 10~ 7 
eV and roughly constant with mass number, while typical RMS matrix elements were on 
the scale of 1 meV. 

Theoretical treatments have been developed that led to model predictions of the 
RMS matrix element based on an effective weak interaction with dominant isovector 
pion and isoscalar p exchange contributions taken from DDH predictions. A recent 
study by Tomsovic et al. (j9lh outlines a statistical spectroscopy approach for interpreting 
the experimental RMS matrix elements based on PV asymmetry data from 238 U and 
104,105, loe.iosp^ targets to set limits on the weak meson-nucleon coupling constants. To 
date, a variety of theoretical approaches have been used to predict RMS matrix elements 
that are in qualitative agreement with experimental data, based on DDH predictions of 
the coupling constants, and which seem to favour a small value of h\ consistent with the 
18 F measurements shown in Figure 121 but a comprehensive program to attempt to analyze 
all the data and set limits on the weak coupling constants has not yet been attempted. 
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2.3 The End of the Road 



Figures El and E| summarize the present state of our knowledge using the meson-exchange 
framework. We have at present reasonably clear experimental constraints on four linear 
combinations of the six weak meson-nucleon couplings from a mixture of few nucleon 
and finite nuclear experiments, not all of which are in agreement with each other. It is 
fair to say that the phenomenological implications that one can draw from these results 
are unclear at best. On the one hand, the pp experiments now yield a rather stringent 
set of constraints on the combinations of the DDH couplings h pp and h pp that govern 
the asymmetry over a fairly broad range of energy. As we discuss below, h pp and h pp 
essentially parameterize the contributions to A z from the lowest-order, short-range PV 
potential in the EFT. To the extent that the energy-dependence of A z is dominated by 
that of the strong interaction phase shifts that enter matrix elements of these operators 
as well as the unpolarized cross section, one need not think of h pp and h pp as being specific 
to the DDH meson-exchange framework. 

On the other hand, the nuclear PV experiments are largely sensitive to two different 
combinations of the DDH parameters shown in Figure |2J h™ c = —{h? + 0.7/i°) and 
h v nc = h\ — 0A2h 1 p — 0.18/1* (0). In the past, it has been the conventional practice 
to project the constraints from the A z measurements onto the h^y plane by using the 
DDH theoretical ranges for the h 1 ^ 2 and . Here, however, we choose not to do so 
since we want to minimize the number of theoretical assumptions used in the extraction 
of information from experiment. Instead, we treat the nuclear experiments separately 
from the pp measurements. In the ideal situation, the analysis of the nuclear experiments 
would yield a self-consistent region for the h'gy - a situation that clearly does not emerge 
from Figure 121 The primary problem seems to be the inclusion of the 133 Cs anapole 
moment constraint, which finds no region of simultaneous consistency with all the other 
nuclear PV experiments. To explain this discrepancy, one might naturally re-examine the 
shell model calculation leading to the cesium band. However, theoretical considerations 
suggest that a more realistic calculation would lead to an even larg er discrepancy 
while the results of more naive shell model computations (BE! 00, HE 
lead to a similar result. Evidently, additional insight into the many-body physics of 
nuclear PV is needed before a consistent phenomenology can be obtained with the meson- 
exchange framework. 

Even if such a consistent picture had emerged from experiment, extraction of funda- 
mental information on the AS* = HWI would still be problematic. To explain why, we 
consider the physics embodied by the hgy. As compared with h pp and hff, these effective 
nuclear couplings correspond to different combinations of the short range EFT operators 
and long-range 7r-exchange PV potential than those that enter the pp asymmetry. How- 
ever, the nuclear matrix elements of these operators sample the spatial-dependence of 
both the operators as well as their action on the nuclear wavefunctions, and we have no 
simple way of disentangling the two as we do for the pp asymmetries. 

To illustrate, we first consider the momentum space form of the p- and tu-exchange 
operators appearing in V^ H . Each term contains a pseudoscalar of the general form 
Bi ■ pj - where the subscripts refer to nucleons i or j - times a function of the momentum 
transfer q = \q\ generated by the meson propagator in the static limit (go = 0) and a 
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hadronic form factor F Pi w (q 2 ) arising from the meson- nucleon vertices: 

^-exchange ~ <?i ' Pj F V tf) [1 + tf/mfa- 1 ■ (7) 

For q << my we may expand in powers of q 2 /m v : 

^-exchange ~ <?i ' P, F V {0) ( 1 + WyF*^) - 1 



in 



v 



where the + ■ • • indicate higher order terms. In a model-independent approach, the 
coefficient [m v F'(0) — 1] of the q 2 /my term in Eq. (jHJ) would be replaced by an a priori 
unknown coefficient whose value would have to be taken from either experiment or a QCD 
computation. A similar statement holds for the higher order terms. At each order n in 
the expansion, one could also include additional pseudoscalar operators having a different 
structure than &i ■pj(q 2 /m 2 / ) n . In effect, the meson-exchange framework imposes model- 
dependent relations between all of the higher-order operator coefficients and those of the 
lowest order terms - relations that may or may not hold in the SM. 

The impact of the higher-order operators on nuclear matrix elements depends on 
both the values of the operator coefficients as well as the spatial dependence of the nuclear 
wavef unctions. Extraction of the hfgy from experimental observables relies on both the 
relationships between these operators implicitly assumed by the meson-exchange model 
as well as on nuclear model-space truncations and other nuclear structure inputs that 
affect the wavefunctions employed. At present, we have no rigorous way to disentangle 
the impact of either on the extracted PV couplings, and the constraints in Figure 121 
may reflect both artifacts of nuclear structure calculations as well as assumed operator 
relations. If one seeks to study the fundamental AS = HWI experimentally, then 
one would like to avoid such an implicit reliance on model-dependent assumptions and 
nuclear structure inputs. The theoretically cleanest way to do so is to exploit effective 
field theory - wherein operator relations are determined systematically from experiment 
- and by studying hadronic PV in few-body systems, for which ab initio theoretical 
computations are available. 



3 Effective Field Theory Framework 

Effective field theories are ideally suited to situations where there exists a distinct hi- 
erarchy of scales. In the presence instance, several scales are relevant: the weak scale 
v = (\/2GfY^ 2 = 246 GeV; the hadronic scale Ahad ~ 1 GeV; the pion mass and decay 
constant, m n w 140 MeV and F n = 93.2 MeV, respectively; and the typical momentum 
Q relevant to a parity-violating hadronic or nuclear process. The distance at which the 
repulsive core of the strong NN potential becomes dominant, r ^ 0.4 fm corresponds to 
a mass scale ^ 500 MeV which, for our purposes, we take to be of order A H ad- The size 
of hadronic matrix elements relevant to the HWI is governed by the the ratio of F 2 to 
v 2 that is typically normalized to the quantity g n as in Eq. 

The remaining scales can be used to construct an effective Lagrangian out of nucleon 
and pion fields, where the operators are organized according to powers of Q/A. For 
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processes in which Q « m^, one should take A = m w , treating the pions as heavy and 
"integrating them out" of the effective theory. For Q ^ m^, the pion must be kept as an 
explicit degree of freedom, and one should take A = Ahad- Consequently, we consider 
two versions of the EFT corresponding to these two different regimes for Q. 



3.1 The Pionless Effective Field Theory 



In the PV EFT without pions, the lowest-order pseudoscalar operators contain four 
nucleon fields and are 0(Q), since they must transform as ff{ ■ pj. At this order, there 
nominally exist ten different contact operators. As shown in Reference (J^J), the most 
general short-range ( "SR" ) potential has the coordinate-space form 



vl 



PV 



SRA 



2 

A3 



Ci + (C 2 + C 4 ) 



r\ + r 2 



C3T1 ■ T 2 +l ab C 5 T?T2 



a x - a 2 ) ■ {-iV,f m {r)} 
T\ + t 2 



+ C! + (C 2 

i (ct! x a 2 ] 

+ (c 2 - c 4 ) 

+C 6 te ab3 



-c 4 
Hv, f w (n) 

Ti - r 2 



C3T1 ■ t 2 +X a6 C 5 rfr 2 



T 1 T 2 l "! + ^2 J 



o\ + ^2) • {-«V, f m (r)} 



(9) 



where T a b = diag(l, 1, 2) and where the subscript "1" on V^ p ^ R (r) essentially indicates 
that this potential appears at 0(Q 1 ) in the EFT. In arriving at Eq. we have intro- 
duced the function f m (r) that is strongly peaked about r = with some width ~ 1/m 
and goes to 5^ 3 '(r) in the zero-width (m — > 00) limit. For practical purposes, we will 
take 0.4 fm. 

At first glance, one sees a dependence on ten a priori unknown constants Ci_6 
and C*i_5 that encode information about the short-distance weak interaction between 
two nucleons 4 . When considering processes with Q « m n , however, not all of the 
operators in Eq. are independent. In this regime, PV observables are dominated by 
mixing between S- and P-waves, for which there exist only five independent spin-isospin 
amplitudes: 



i) d t (k), representing 3 «Si(/ = 0) 

ii) d®' 1,2 (k), representing 1 S (I = 
operators respectively; and 



- 1 P±(I = 0) mixing 

1) — 3 Pq(I = 1) mixing generated by / 



0,1,2 



111 



c t (k), representing 3 Si(I = 0) — 3 Pi(I = 1) mixing 

where we have used the notation of Danilov (|92l ) and Desplanques and Missimer ffiol 
liil ). At the low energies relevant to the pionless EFT, the energy dependence of these 

The combination C2 — C4 does not appear, so that only ten combinations of the eleven constants 
Ci_6 and Ci_5 appear in Eq. ©. 
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amplitudes is dominated by the strong interaction phase shifts. Denoting the spin singlet 
and spin triplet strong interaction S-wave scattering amplitudes as m s (k) and m t (k), 
respectively, we have 



d t {k) 
c t (k) 



X t m t {k) + ■ 
Km s {k) + ■ 
p t m t (k) + • • 



(10) 



where the + • • ■ indicate small corrections to the energy dependence arising from the 
strong P-wave phase shifts. In his early formulation of the problem, Danilov (|92l ) ar- 
gued for the forms in Eq. (fTU|) . omitting the small corrections indicated. In the EFT 
framework, it is straightforward to derive the proportionality of the S-P amplitudes and 
the rrii(k) by computing the relevant T- matrix elements and summing up the strong 
rescattering contributions. In doing so, one finds that at low energies where P-wave 
rescattering contributions are small, the Aj and pt are given by the ratio of the lowest 
order S-P scattering amplitude to the lowest order, parity-conserving S-wave amplitude 
in a given spin-isospin channel (@). 

The coefficients A(, A*, and pt are themselves proportional to various combinations 



of the Ci and Ci appearing in V^gp^f). In the zero-range 

A* oc (d - 3C 3 ) - {Ci - 3C 3 ) 
X° s oc (d + C 3 ) + (Ci + C 3 ) 
X] oc (C 2 + C 4 ) + (C 2 + C 4 ) 

A? oc • 



m 



oo) limit, one has 5 



:ir: 



}/3(C B + C B ) 



1 



Pt oc -(C2-C4) +C 6 . 

In effect, at low energies, five of the operators in Eq. © become redundant, leaving only 
five independent S-P amplitudes. 

Inclusion of finite-range effects leads to modifications of these relations. Arriving 
at exact relations using state-of-the-art NN potentials remains an unfinished task for 
many-body theorists. However, we provide approximate expressions by drawing on the 
work of Desplanques and Benayoun (joih . We obtain 



m N p t = 


2 


B 2 


m N X t = 


2 


~B 4 


m N X° s = 


2 


B 6 


m N Xl = 


2 


B 6 


m N X 2 s = 


aVq 

A 3 


B 6 



~:Co 



;C 4 +C 6 )+B. 



Co 



~Ca — Cf 



(12) 



5 The last equation of Eqs. (46) in Reference (|3) contains an error. The sign in front of C 6 should be 
"+" rather than "-" . 
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where A = m^w? p j 'A 3 and where the Bk are linear combinations of the (3fj of Refer- 
ence (jjih. For example, using the values of those constants obtained with the Reid Soft 
Core (RSC) potential, we obtain B k = (-0.0043, 0.0005, -0.0009, -0.0022, -0.0067, 0.0003) 
for k — 2, . . . , 7, respectively 6 . 

The coefficients -82,5,6 multiply the combinations of Cf and Ci that one expects to 
arise in the lowest order EFT [ Eq. The remaining terms are generated by finite- 

range contributions that occur beyond leading order. Although we have not included the 
full set of PV operators and amplitudes that occur at next-to-leading order (NLO), the 
magnitudes of the -83,4,7 - relative to the -82,5,6 - give an indication of the magnitude of 
higher-order effects and of the error associated with working to lowest order. The RSC 
values obtained in Reference f)94T ) give \B 3 /B 2 \ = 0.12, \B 4 /B 5 \ = 0.41, and \B 7 /B 6 \ = 
0.04, suggesting that the impact of neglected higher-order contributions are generally 
small except in the case of X t . 

Before considering the application of this framework to specific observables, it is 
useful to obtain theoretical predictions for the quantities p t and A Sjt . To that end, we 
first consider the correspondence with the DDH meson-exchange model and delineate the 
relationship between the Ci and Ci and the DDH parameters. Using 



fm(r) 



m 
Arcr 



exp(— mr) 



(13) 



with m being the parameter that defines the range of the PV potential, and letting 

A 3 



A 3 



m N mf M 



(14) 



for M = p, u we have 7 



C{ 



DDH 



qDDH 



c 



DDH 



1 

~2 

1 — 



qDDH 
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riDDH 
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iA 3 h 1 

2 UjiJ^ Li) 



- ! A 3 o h 1 



(15) 



and 



C 



DDH 



(JDDH 



C 



DDH 



c 



DDH 



1 + Xu 



1 + 



1,2, 



3-5 



(16) 



6 We emphasize that the expressions in Eq. (|12J) are applicable to the EFT without pions. For a 
discussion of the modifications due to inclusion of explicit pions, see Sec. 13.21 below. In particular, the 
lowest-order contribution from single pion exchange appears in p t . 

7 The corresponding expressions given in Eq. (141) of Reference (@) contain typographical errors. The 
quantity Ay should contain only two powers of tum in the denominator and each of the Ci should be 
proportional to the product of a strong coupling gM and the relevant h l M as in Eqs. (|14I16(I here. 
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Using Eqs. ()15I16J) one can easily obtain the expressions for p t , X t , and A°' 1,2 and employ 
the DDH best values and reasonable ranges for the PV meson-nucleon couplings to obtain 
the predictions for the PV LECs listed in Table |2j The results in columns 2-4 were 
obtained by retaining only the combinations of the Ci, C{ that arise at lowest order in 
the EFT. Column five contains estimates of the size of higher-order contributions, based 
on the DDH best values for the Cj, Ci and the the terms in Eqs. ()12|) proportional to 

-53,4,7- 

To obtain a sense of the possible variations in theoretical predictions for the PV 
LECs from their correspondence with the DDH parameterization, we consider two ap- 
proaches. First, we vary the values of the strong couplings in Eqs. ( 115)1 in accordance 
with the Bonn one pion exchange potential as suggested by Miller (J95J). As indicated 
in Table 121 doing so changes both the overall magnitude of the Cj and Ci as well as 
the relation between the two and leads to generally wider ranges for the PV LECs than 
obtained with the values of the strong couplings originally used by DDH. Second, we 
give expectations using naive dimensional analysis (NDA) considerations, as discussed in 
Section 13.21 below. The NDA arguments suggest that the magnitudes of the Ci and (7, 
ought to be of order 167r 2 ~ 150, but do not fix the signs of the Cj and Cj. To translate 
the NDA estimates into predictions for the five LECs, we simply take the magnitudes of 
the combinations of Ci and Ci appearing in Eqs. f[T2|) to be the NDA expectation for any 
one of them. 



Table 2: Predictions for the five, PV low energy constants (LECs ) characterizing hadronic 
PV in the pionless EFT. All values are quoted in units of q n = 3.8 x 10~ 8 . Estimates 
are obtained using Reid Soft Core potential as in Reference ; DDH values are taken 
from Reference \Q) 



PV LEC DDH Best DDH range DDH plus Bonn higher order NDA 

m N p t 0l)5 0.07 -»• 0.03 0.18 -»• 0.08 ±0.006 ±0.38 

m N \ t 0.84 -1.00 -> 2.48 -2.44 -»■ 6.12 ±1.2 ±0.31 

m n \° s 3.82 -4.86 -> 11.7 -9.84 -> 22.96 ±0.09 ±0.64 

m N Xl 0.37 0.64 -> 0.21 1.53 -> 0.54 ±0.0006 ±0.64 

m N \l 2.72 2.17 -> 3.15 3.83 -> 5.55 ±0.06 ±3 



We note that in obtaining the correspondence between the PV LECs and the DDH 
predictions in the meson-exchange model, we have not included contributions from the 
parameter h 1 ' that appears in V^ u . Using the estimate of Reference for this param- 
eter would increase the magnitude of p t appearing in Table El though not substantially. 
Generally, analyses of hadronic PV using the meson-exchange model have neglected this 
parameter since its contribution to has the same spin-isospin structure as for the 

7r exchange contribution but is suppressed by its short range. In the EFT framework, 
this term corresponds to the operator proportional to Cq that - along with C*2,4 - con- 
tributes to p t . Since the coefficients of 6*2,4,6 in Pt have comparable magnitude, we see no 
model-independent reason to neglect the Cq contribution in the most general analysis. 

With the foregoing set of benchmarks in hand, it is instructive to consider the 
dependence of various few-body PV observables on the \ and p t and to outline a program 
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of high precision measurements that could be used to determine these parameters. The 
few-body PV observables of interest include: 

• Polarized pp scattering at 13.6 and 45 MeV, yielding the asymmetry A pp , 

• Polarized pa scattering at 46 MeV, giving A pa , 

• Radiative np capture at low energy: np —>■ d'j, yielding the photon asymmetry A d 

• Radiative np capture with unpolarized neutrons, giving the photon circular polar- 
ization P d or alternatively, the asymmetry A\ in 7<i — > np, 

• Rotation through an angle <fi about the momentum direction of polarized neutron 
spin passing through 4 He, from which one extracts the quantity d<j) na /dz, 

• Radiative capture of polarized neutrons on deuterium at threshold, nd — > tj, yield- 
ing the photon asymmetry A t . 

Explicit expressions for these quantities in terms of the S-P amplitude parameters 
have been given in references (jiol I93T ) and elsewhere. In Table El we give the coefficients 
of the five PV LECs as they appear in various observables. Theoretical expectations 
for these observables in the pionless EFT can be obtained using Tables 121 and EJ We 
emphasize that these expectations will, in general, differ when going to the EFT with 
explicit pions discussed in Section T3.2I In particular, the parameter m^Pt that governs 
the asymmetry A d will be dominated by LO pion exchange, assuming h\ has its natural 
size. In Table we also show the most precise experimental results that have been 
published to date for these observables. Notably, only A pp and A p z a have been measured 
to sufficient precision to establish a nonzero PV effect that can be used to constrain 
the PV LECs. A review of recent, ongoing and prospective efforts to obtain precision 
measurements of these important few-body PV effects appears in Section l3~4l below. 8 

It should be remarked that the l ong itudinal asymmetry A p z d in pd scattering was 
measured to high precision at 43 MeV (|98l. IQ9T) in the 1980's, but remains to be analyzed 
in a theoretical framework accounting for PV in both elastic scattering and breakup 
channels, both of which contributed to the experimental signal. There exist also several 
additional possibilities for few-body experiments, for which we have not yet obtained ex- 
pressions for the PV observables in terms of the PV LECs. These include measurements 
of PV neutron spin rotation on hydrogen, d<p np /dz and deuterium, d<p nd /dz, and one 
could also consider the circular polarization P* in nd — > £7 as well as perhaps even the 
transmission asymmetry of unpolarized neutrons through polarized 3 He. With realistic 
prospects for performing some or all of these measurements in the future, deriving the 
appropriate expressions would clearly be important. 



8 The measurement of is particularly challenging because of the limited sensitivity of conventional 
7 ray circular polarimeters in the ~ MeV energy range; an alternative measurement of JC L close to 
threshold yields the same physics and may be more accessible to experiment. A possible experiment 
is at the early stages of development in Athens fTf!) . as mentioned in Section WJA but is not discussed 
further in Section ITU 
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Table 3: Sensitivities of selected PV observables to the five PV LECs. The first column 
gives the observable, while subsequent columns give the coefficients of a given PV LEC. 
For the pp asymmetry, k is the incident proton momentum in the lab frame. The final 
columns give the most precise experimental limits to date and their references; for the pp 
case, we quote the 13.6 and 45 MeV measurements, to which the lowest order EFT best 
applies. Note that the pa asymmetry is evaluated at 46 MeV; d(p na /dz is evaluated in 
rad/m. 
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Limits of the pionless EFT 

Before discussing the EFT with pions, it is instructive to investigate the limits of 
validity of the pionless theory by considering the the low-energy pp asymmetry as an 
illustrative example. Again, considering only S-P mixing, it is straightforward to show 
that @): 

^^jtReKy)]^^ (17) 



<T4- + CT- 



\m s (k)\< 



where where we have neglected small corrections arising from the P-wave phase shifts 
as before. A more complete decomposition of AT in terms of higher partial waves was 
first worked out by Simonius ([461 ) and subsequently studied by several authors. A recent 
analysis using state-of-the-art NN potentials was performed in Reference (jH3) and used 
to extract the DDH parameters hP^ and hT ■ 

Following Simonius' original formulation, one may write the asymmetry as 

AT= E fME)KME,9) , (18) 

even J 

where the "±" indicate the orbital angular momentum L — J ± 1 and where E and 9 
are the energy and scattering angle, respectively. The fj± are "reduced" PV transition 
amplitudes and the Kj± contain all the dependence on the strong phases that arise from 
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rescattering. For the J = partial wave, only / + exists, and it is proportional to the 
combination d pp of the S-P amplitudes d\ that enter the pp process: 



df = -2z/ 0+ exp i [6(%) + 5( 3 P )] . (19) 

To the extent that one may neglect the finite (but short) range of the PV potential, d pp 
and / + are proportional to \ pp . Moreover, in obtaining Eq. (|17|) . we have neglected the 
dependence of d pp on <5( 3 P ). Doing so is equivalent to taking K(E, 6) oc k. 

While this approximation holds to a high degree for accuracy for low energy inter- 
actions, it breaks down for E ~ 100 MeV (k ~ 300 MeV). Due to cancellations between 
the effects of the S- and P-wave phase shifts in d pp that occur above this energy, the con- 
tribution of the J = partial wave to A pp falls rapidly, going to zero at E — 227 MeV. 
Moreover, the contribution from the J = 2 partial wave - dominated by the 1 D 2 - 3 P2 
(/2-) mixing - becomes appreciable. Taking advantage of the former, the beam energy 
for the TRIUMF A^ p (221 MeV) measurement was optimized, accounting for finite accep- 
tance of the apparatus, to ensure that that the /o+ contribution was entirely cancelled, 



so that the asymmetry was determined almost entirely by the j^- contribution ([4 71 154I). 
Thus, a combined analysis of the TRIUMF and lower-energy A pp measurements yields 
constraints on the two amplitudes fo+ and /*2-. In effect, one may treat the values of h pp 



and h pp obtained in Reference (54) as equivalent parameterizations of these two partial 
wave transition amplitudes and need not tie them specifically to the meson-exchange 
framework. 

Clearly, at the energies of the recent TRIUMF experiment, the lowest-order EFT 
is no longer applicable, and one most include 0(Q 3 ) operators in that characterize 
P-D mixing. Doing so introduces a host of new, a priori unknown operator coefficients. 
At the same time, the redundancy of operators in the 0(Q) in Eq. (J2J) that holds for 
the S-P amplitudes breaks down, and all ten lowest order operators become independent. 
Determining all of these constants from experiment would be unrealistic, so we will 
restrict our attention to the energy range where the lowest order EFT applies. 

3.2 PV EFT with Pions 

For PV processes involving few-body nuclei, the relevant energy scale Q is no longer set 
solely by experimental kinematics, but also includes the relevant internal momentum of 
the bound nucleons. Since the latter can be as large as the Fermi momentum of ^ 200 
MeV, it is no longer reasonable to treat the pion as heavy. For decades, phenomenolog- 
ical strong interaction potentials for light nuclei have included a long-range 7r-exchange 
component. In the case of the EFT formulation, treating pionic contributions consis- 
tently has presented challenges. The difficulty arises from the presence of two-nucleon 
poles in iterated 7r-exchange amplitudes whose contributions are enhanced by ~ m^jQ 
relative to naive expectations. These enhanced contributions spoil the "power counting" 
in Q/A that is essential to the success of the EFT approach. Thus, one must sum strong 
pion exchange to a ll orders in order to obtain a consistent treatment. Indeed, as shown 
in Reference treating the pion perturbatively as in the framework of references 



fjlOlL 11021 ) does not lead to a convergent expansion in all channels of the NN interaction. 
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An alternate formulation, originally proposed by Weinberg 

Jiollioi), entails per- 



forming the all-orders resummation in terms of the effective strong potential. In order 
to be self-consistent, however, it appears that one must simultan eously perform the EFT 
expansion of the potential about the chiral limit ( 10,4 lOfil . llOTf ) : m w = 0, since a full re- 



summation of the chiral symmetry-breaking co mpon e nt o f the one pion-e xchange (OPE) 



potential leads to inconsistent renormalization (101, 102. Il05l llOo . Il07f ). Although the 
chiral expansion treatment of the Weinberg approach is still under development, we follow 
Reference and employ it here. 

The basis for the EFT with pions is chiral perturbation theory (%PT) whose formal- 
ism is well-known and will not be repeated here. However, we note that its form ulation 
for hadronic PV processes was first written down by Kaplan and Savage (|l08| ). whose 
notation follow. The non-linear dependence of the effective Lagrangian on the pion is 
implemented via the field 



where ir a , a = 1, 2, 3 are the isospin components of the pion field. Defining the quantity 

X 3 =£t r 3£_£ T 3 e t ; (21) 

the lowest-order Lagrangian for a PV interaction of the pion with a single nucleon is 



= -7fif x - N 

= —ih\{pn-K + — np-K~) + • • • (22) 

where the "+ ■ • ■" indicate the higher order terms in odd powers of (7r a T a / F n ) that arise 
from expanding the exponential in £. The leading term in the Lagrangian in Eq. ()22j) is 
identical to the PV Yukawa interaction in the DDH model. Consequently, its contribution 
to the PV NN potential will be the same as the first term in Eq. (J2J): 



t/PV ( ^_ j h l9Am N 

V {-1, LR)V> — 1 




Pl~P2 f s 

— ,w v {r) 

ZrriN 



(23) 



The "—1" subscript indicates that this long-range ("LR") potential is 0(Q 1 ), a feature 
most readily seen from its momentum-space form: 



v pv _ _■ 9aK ( ri xr 2 \ (a, + a 2 ) ■ q 



where q = pi — f>{ = p 2 ' — p 2 is the three-momentum of the exchanged pion, p^ (p/) is 
the initial (final) momentum of nucleon i, and q = \q\- 

Since the PV potential must transform as a pseudoscalar, the operators in it will 
contain odd numbers of derivatives. Thus, one would expect the sub-leading terms to 
be 0(Q), as in V^ V SR ). In principle, loop corrections to the Lagrangian (J2~2~j) or to 
VZY LR)(r) could bring in a factor of m n or p 2 /m n , leading to an O(Q ) component. 
Explicit computations, however, indicate that no such contributions exist. Consequently, 
the subleading components of the potential start off at 0(Q) (next-to-next-to-leading 
order, NNLO), and it is convenient to distinguish them according to their range: 
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i) Short range: V^s R) as in Eq. Q but with A = A H ad ~ A x = AttF^. 

ii) Medium range: V^j* mr n, generated by the two pion-exchange diagrams of Figure 
d?) and proportional to h\. The structure of the operator is most conveniently given 
in momentum space, as it carries a non-analytic dependence on pion momentum 
and mass 



where 



and 



^(i?mr)(£) 



L{q) 
H(q) 



t? + r 9 2 



-i (cri x a 2 ) -q 



+C 2 "(q)ie ab3 [f 1 x r 2 ] 3 {B x + a 2 ) ■ q) , 

Cf( ff ) = AV2irg\hlL{q) 
-V2ng A hlL(q) + ^tt [3L(q) - H{q)] g\h\, 



4m 2 + q 2 

q 

4m 2 



In 



' 4m 2 + g 2 + g 
2m„- 



4m 2 + g 2 



L(g). 



(25) 



(26) 



(27) 



iii) Long range: ^ V LR ), generated by one- loop corrections to the PV itNN Yukawa 
and parity conserving strong vertices as well as by additional operators having a 
distinct structure from the lowest order potential. As discussed in Reference (0), the 
impact of all but one of the NNLO PV operators can be absorbed in to VF^_ LR \ and 
^Ol V sr) through a suitable redefintion of the operator coefficients. The remaining 
PV operator and PC operators give rise to the momentum space potential 



VTlr(Pi 



P 2 



g A K a ( f x x f 2 



A F 2 



(Ti-Pi x p x a 2 ■ qi 
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9aK 
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(|pi| 2 + |p 1 / | 2 )a 1 -g+(1^2) 
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+4 [<?i ■ Pi q- pi + Bi- Pi q ■ pi + (i 



2)] 



(2f 



where qi = p/ — pi and k\ a is a constant that must be determined from experiment. 

In addition to considering the potential through 0(Q), one must also include two- 
body current operators that contribute to the same order when considering PV processes 



36 



involving photons. These operators include the standard PV meson-exchange currents 
associated with the 7r-exchange potential and those arising from the covariant derivatives 



J{x u x 2 , q) = V ZXZ^ e*** t+t 2 a.xqa.-r H^r) + (1 <-> 2), (29) 



and CV is an additional LEC parameterizing the leading PV NNif-f interaction. 

PV EFT with pions: new features 

In comparison with the DDH meson-exchange potential, the EFT with explicit pions 
introduces several qualitatively distinct features. First, the 0(Q), long range operators 
in Eq. (|28I29J) have no analog in the DDH framework and introduce two new unknown 
constants, k\ a and C % . Given the novel nature of these operators, the impact of their 
contributions to the observables discussed above has yet to be determined with explicit, 
few-body computations. Based on the isospin structure of these two operators, however, 
one would expect them to contribute to the same processes, such as np — > dj, that are 
sensitive to the LO 7r-exchange potential. Power counting implies that the magnitude of 
their contribution should be smaller than that of V^f LR ), but their long range character 
suggests that they should have a greater impact than the operators in V^ V SR ). Future 
few-body calculations should test these expectations. 

The two 7r-exchange (TPE) medium range potential, V^ V MR) , is similarly a new 
feature of the EFT framework. In the past, oth ers h ave attempted to introduce PV TPE 
using model frameworks (see, e.g., Reference (jl09j )). but to our knowledge, the result 
in Eq. (J25)) gives the first formulation that is model-independent and consistent with 
the symmetries of QCD. The operator coefficients Cf 71 and Cg 71 " have been labelled to 
indicate their correspondence with the operators in V^gp^, but as indicated in Eq. (J2HJ), 
the coefficients are fixed in terms of h\ and are not independent free parameters. The 
operator proportional to Cf^ has the same isospin structure as V£Y LR ) and will contribute 
to any process - such as np — ► dj - that is sensitive to the LO 7r-exchange potential. The 
operator proportional to has the same structure as the /z^ (1 + x<J) term in the DDH 
potential and will, therefore, generate a medium range contribution to any observable 
sensitive to the latter combination in the meson-exchange model. In particular, this 
operator will contribute to both the f + and / 2 - partial wave terms in A p z p , implying that 
this observable is sensitive to h\ at the same order in Q as the the short-range effects 
enter. In either case, proper inclusion of TPE medium range potential will affect the 
determination of h\ obtained from a global analysis of few-body PV observables. 

In principle, one should also take into account three-body PV forces that arise in 
the EFT with pions, as we will discuss few-bodv experiments in systems involving three 
or more nucleons. As discussed in Reference (|8j), three-body PV forces do not arise at 
0(Q), so we restrict our attention to the two-body PV EFT interaction. 




where 




(30) 
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PV LECs: naive dimensional analysis 

Beyond these qualitative observations, we are not able to make any quantitative 
statements regarding the relative importance of the new features of the PV EFT with 
pions. Evaluating their contributions to specific observables now constitutes an open 
problem for few-body theorists. That being said, it is instructive to estimate the size 
of the constants Cj, Ci, h\, k\ a , and C n that arise in the EFT. A systematic way of 
doing so - kn own as "naive dimensional analysis" (NDA) - was developed by Georgi 
and Manohar and has successfully explained the size of a variety of LECs in x?T. 

According to NDA, one should construct operators by scaling fields and derivatives to 
their natural scales: 

/ n \ 1 / „ \ V / AT AT \ f/ 2 

x (A x ig 2 x (g n ) n , (31) 

where d, p, f = 2k, k and n are positive integers and where g n is as given in Eq. (jlj). 
Consequently, one expects the PV LECs to have the magnitudes 

hi ~ (^)^~10^ (32) 
A/ 2 




Cud ~ g v ~\Wg v (33) 

kl\ a = g« ■ (34) 

It is interesting that since A x /F n = 4tt ~ 12, the NDA estimates for hi and the Ci, 
Ci are roughly equal to the expectations based on correspondence with the DDH best 
values. In contrast, one expects the new long range LECs k\ a and C n to be an order 
of magnitude smaller than h\. This difference, however, simply reflects the conventions 
used above in normalizing the various operators. 

PV EFT: phenomenology 

The phenomenology of the PV EFT with pions is clearly more challenging than for 
the pionless theory, since one encounters additional unknown parameters when working 
to 0(Q), and since a variety of new contributions remain to be computed. In principle, 
there exists a viable experimental program that could determine the constants to this 
order, including the six measurements highlighted in Section I3~T1 for the pionless theory 
plus two additional, independent experiments. As noted in Section I3~T1 possibilities for 
the latter include a careful analysis of the pd scattering asymmetry results for A p z d (98|, 9i 



and future neutron spin rotation measurements on hydrogen or deuterium. Recently, the 
possibility of performing PV photo- and electro-production experiments on the single 
nucleon to determine the of ttN and 7iV couplings has received considerable attention. 
Carrying out these experiments - which we discuss below - would provide additional, 
independent input for the determination of the PV LECs. 

Theoretically, the extraction of these constants from experiment will require new 
calculations to determine the contributions from (a) the medium-range, TPE potential 
(J2HJ), (b) the NNLO single pion-exchange potential (|28|l. and (c) the meson-exchange 
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current ([2 9ft. At present, only the dependence on h\ generated by LO pion-exchange is 
known: 

m N p t = 1.04/4 + m N p* R + m N Ap t , (35) 

where gives the dependence of the 3 Si- 3 Pi mixing on the constants appearing in 
Vi sr as m Eq. JE}; Ap t gives the presently unknown contributions generated by Vimr 
and Vi\ K ; and the first term on the right side of Eq. (J33j) is generated by LR as 
computed by Desplanques and Benayoun using the Reid Soft Core potential (jjih . Since 
^i P Xir contains spin-isospin structures corresponding to both the C% and C2 terms in 
Vig R , the LEC \] will also contain a dependence on h\ generated by two-pion exchange. 
In addition to computing these new contributions to p t and Aj, theorists must also de- 
termine new meson-exchange current contributions to processes such as n + p — > d + 7, 
generated associated with V['\ lR and V±\ R and required by gauge invariance, as well as 
the contribution from the new current in Eq. (|29j) . 

Assuming that a successful program is completed and the complete set of PV LECs 
through 0(Q) are extracted from experiment, the values of these parameters would the 
provide model-independent benchmarks for Standard Model theory. In this case, the 
theoretical challenge would be analogous to the one encountered with ;\PT for pseu- 
doscalar mesons, where for example, at 0(Q 4 ), there exist ten independent LECs that 
have been determined from experiment. The theoretical task is now to explain how the 
dynamics of QCD give rise to the values of these constants, and to that end, a number 
of approaches have been pursued. Ultimately, of course, one would like to compute these 
constants using lattice QCD, but given the difficulties in putting two or more hadrons 
on the lattice, approaches based on symmetry arguments or models are an attractive, in- 
terim alternative. A particularly fruitful direction involves taking the limit of QCD with 
a large number of colors (Nc), wherein one expects the exchange of heavy mesons, such 
as the p and lu, to dominate the underlying QCD dynamics of the LECs. In the Nc = 3 
world that we inhabit, this large- Nc picture of "resonance saturation" works remarkably 
well in accounting for the values of the constants. It remains to be understood why large 
the large Nc limit is so successful in this case, and future lattice QCD computations 
should address this problem 9 . 

In the case of the PV LECs, the results from experiment should teach us whether 
the large Nq resonance saturation picture applies to the AS* = HWI involving baryons 
as well as to strong interactions between light mesons. In effect, the DDH model as- 
sumes the validity of resonance saturation, albeit with a truncated spectrum of exchanged 
mesons that may or may not reflect accurately the underlying dynamics. As discussed 
at the outset of this article, there exists ample evidence that QCD symmetry arguments 
fall short when confronting the phenomenology of the AS* = 1 HWI, so one should apply 
caution when adopting another one (viz, large Nc) to predict weak, hadronic AS = 
processes. At the same time, one would like to derive as much model-independent infor- 
mation as possible on the HWI in each sector, so that one can gain new insights into the 

9 At the same time, it has become important to know the values of certain 0(Q 6 ) constants that 
presently cannot be taken from experiment and that are needed for the extraction of the Cabibbo- 
Kobayashi-Maskawa matrix element V us from K e ^ decay data. The insights and techniques developed 
to explain the 0(Q 4 ) constants will be essential in obtaining reliable theoretical values for the unknown 
higher order terms. 
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puzzles associated with strangeness-changing processes. For example, if one ultimately 
found a set of PV LECs that agreed with expectations based on NDA and resonance 
saturation, one might conclude that the breakdown of symmetry-based expectations in 
the AS = 1 sector is associated with the dynamics of the participating strange quark. 
On the other hand, should the PV LECs depart substantially from NDA and large Nc 
expectations, one would look elsewhere to determine dynamics general to all sectors of 
the HWI. 



3.3 Recent Theoretical Work 

As the foregoing discussion makes evident, 
theoretical work within the context of the 
seen initial efforts in this direction, and we 



there now exists ample motivation for new 
EFT for hadronic PV. The past decade has 
review some of this work here. 



Few-Body Systems. 

In the two-body sector, recent interest has focused on the asymmetry for np — > 
g?7, where new computations using EFT and Green's function Monte C arlo m ethods have 
been used. The lowest order EFT computation yields the asymmetry (jlllh 

2m N Re[(X + Y)*W] 



r 



2\X\ 2 + \Y\ 



(36) 



where X and Y give contributions to the parity-conserving np — > d'y amplitude for an 
initial 3 Si and 1 S state, respectively, and W gives the PV amplitude: 
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(37) 



with 7 = y/m^B and B being the deuteron binding energy 10 . From these expressions, 
one obtains = — 0.17/i*. The coefficient of h\ in this result is nearly a factor of two 
larger than in previous, wavefunction-ba sed compu tations, and stimulated considerable 
follo w-up theoretical activity (J41 

, H M M m In Particular, the authors of Refer- 
ence (|ll4f ) computed the asymmetry to NLO in the Weinberg scheme, using two-body 
wavefunctions derived from the Argonne v±s potential, and obtained A^ = — 0.10ft.*, in 
close agreement with previous results obtained using Siegert's theorem, A^ ~ — O.ll/i*. 
Subsequently, the authors of Reference ()4lf ) performed a wavefunction-based computa- 
tion in the DDH framework, using Argonne v is, Nijmegen-I, and Bonn-CD wavefunctions 
and found A* = -(0.106 0.109)/^ + ■■■, where the range corresponds to the choice of 
different potentials and where the + • • ■ indicate small, O(l0~ 9 ) contributions from the 
short-range terms in the DDH potential. 

From a more academic perspective, several computations of the deuteron anapole 
mom ent have been performe d usi ng EFT and wavefunction methods 
118). The LO EFT result is Elf 



T15L ITTgI. (117 



eg A h\m 
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«i- 



m n + 7 
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2m^ + 97 



(3f 



(m^ + 27) 2 6(m^ + 2 7 ) 2 _ 

10 I n writing down Eq. I|37|l . we used an overall sign that is opposite t he one ap pearing in Reference 
(|X X if) . thereby following the conventions used elsewhere in the literature l)ll2llll^) . 



40 



where K\ ~ 1.85 is the isovector anomalous magnetic moment of the deuteron. The nu- 
merical value of obtained from this expression has the same sign but a magnitude 
that is ~ 40 — 50% larger than results obtained using a wavefunction computation (j^sh or 



Weinberg EFT approach (jll7T ). From a practical standpoint, the impact of the deuteron 
anapole moment would be most relevant to the interpretation of PV elastic ed scatter- 
ing, wherein it would generate a potentially important contribution to the isoscalar axial 
vector response (j35h. Since the latter vanishes at tree-level in the Standard Model, it 
is particularly transparent to higher-order effects, such as electroweak radiative correc- 
tions, the strange quark axial vector current, and hadronic PV. To date, however, no 
experiments have been proposed to study the PV elastic deuterium asymmetry. 

Single Nucleon Sector. 

From the standpoint of theoretical interpretability, PV pion photo- and electropro- 
duction processes involving single nucleon targets offer several advantages. In particular, 
they provide a means for accessing the PV nNN couplings directly without having to 
disentangle the short- and medium-range effects discussed above. Moreover, the use of 
%PT to describe low-energy pion-nucleon interactions is well established and has been 
thoroughly studied in the parity conserving sector. Looking to future work in QCD, it is 
likely that attempts to compute the PV 7rNN couplings on the lattice will precede any 
efforts to study the short range PV interaction. 

Historically, single nucleon, PV pion photo- and electroproduction processes were 
studied in t he m eson-exchange framework two decades ago by Woloshyn and Li 

and Henley (|l20[ ). These authors found that one should expect PV asymmetries for the 
scattering of longitudinally polarized photons to be of order a few x 10~ 7 assuming the 
DDH best value for h\, while those for electroproduction could be up to two orders of 
magnitude larger. Given the experimental challenges associated with such tiny asym- 
metries, and the prospect of measuring considerably larger effects in light nuclei, the 
prospects for carrying out single nucleon studies were largely ignored for many years. 

Recently, however, advances in experimental techniques for measuring (9(10~ 7 ) 
photo- and electroproduction asymmetries have stimulated renewed interest in this di- 
rection. Theoretically, Chen and Ji refor mula ted the earlier work for near threshold PV 



pion photo- (jl2lt ) and electroproduction (J122J) using heavy baryon yPT. Sub leading con- 
tributions were subsequently considered by the authors of Reference (|l23| ). To 0(Q), 
one has for the threshold photoproduction asymmetry 



V2F n 
9Am N 



Hp - H n 1 H 



hi + i^pe, (39) 

9aK 



where A x = AnF n and the terms proportional to m n give the 0(Q) contributions. A 
measurement at Jefferson Laboratory could in principle yield a result for 5 7 with sta- 
tistical accuracy at better than the 10~ 7 level; however, substantial technical challenges 
associated with the implied current mode pion detection would have to be overcome to 
design a successful experiment. 

Considerably larger photo- and electroproduction asymmetries may be observed in 
the vicinity of the A (1232) resonance. The asymmetry in this region is dominated by 
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the l owest order PV jNA interaction that does not contribute strongly at lower energies 

AN-y ■ e 



urn- 

*-PV ~~ b ' 



d\ A+ lx p + d A A° 7A n F» x + h.c. (40) 



where the are low energy constants that govern the strength of the PV transition. 
For « rriA — mjv, the PV photoproduction asymmetry is 

B^-*$T + - (41) 



where C3 ~ 2 is the transition magnetic moment and where the "+ • • • " ind icate higher 
order, chiral corrections. The latter have been computed in Reference (jl25|) and shown 
to be relatively small. Thus, the size of the resonance asymmetry is essentially set by 

4- 

The authors of Reference (|125T ) noted that a determination of d% could provide 
additional insights into the puzzles surrounding the AS = 1 HWI discussed earlier. In 
particular, if the dynamics responsible for the enhanced PV hyperon radiative decay 
asymmetries also occur for the AS = PV N —>■ A transition, then one might expect 
as large as a few xlCT 6 . Such an effect could be observed in forward angle PV 
electroproduction experiments, for which the contribution of Z° ex change becomes kine- 



matically suppressed, thereby exposing the g?a contribution (|125f l. On the other hand, 
if the occurrence of enhanced PV asymmetries is unique to the AS = 1 sector and is 
associated with presence of valence strange quarks, then one would expect to smaller 
by an order of magnitude or more. 

With this motivation in mind, possibilities for measuring B^ using both the GO 
and Q W eak instrumentation at Jefferson Lab are being actively pursued. The GO collab- 
oration will measur e inc lusive pion asymmetries in an upcoming backward angle run on 
a deuterium target (|l26| ). A measurement of the PV asymmetry in inclusive inelastic ep 
scattering at much lower Q 2 to ~ 0. 09 p pm is also envisioned as a future enhancement 
of the Q W eak experimental program ()127l ). 

As an alternative to photo- and electroproduction, one may also consider PV Comp- 
ton scattering from the nucleon. Computations of the asymmetry for scattering with 
eithe r pol arized protons or polarized photons have been carried out in references f)f28| ) 
and (|l29T l. The asymmetry in both cases is proportional to /ij at leading order. For 
E 1 << and center of mass scattering angle 6 = n/2 one has (jl29) 

A, p ^ p ~ 8.8 x 10- (-^) (J^j) 3 (42) 



5 x 10- 7 J V70 MeV 

where one expects higher order corrections to yield corrections of order 25% and where 
hi has been scaled to the magnitude expected from NDA. Thus, one could expect to 
see an asymmetry of order a few x 10~ 8 - roughly the size of the asymmetry expected in 
n+p — > d + '-f. The magnitude of the asymmetry for Compton scattering with a polarized 
target is similar. To date, no experimental proposals have been developed for measuring 
either asymmetry. 
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Computing h\ in QCD 

Much of the recent theoretical focus has fallen on extracting the leading PV pion- 
nucleon coupling in a way that does not require knowledge of the other PV LEC's or 
of many-body nuclear physics. The interest in h\ has also stimulated new analyses of 
QCD predictions for this quantity that attempt to go beyond the work of DDH. A first 
principles computation will ultimately require use of lattice QCD, and while we are not 
aware of immediate plans to carry out such a calculation, some of the necessary theoretical 
groundwork has been recently laid. In particular, tractable lattice computations typically 
involve use of quarks that are heavier than the physical light quarks, so in order to obtain 
a physically realistic QCD prediction from a lattice result, one needs to know the quark 
mass dependence of a given quantity 11 . To that end, provides the necessary link, 
since the chiral expansion in powers of is equivalent to an expansion in ^/ m q . 

The first such analysis of h\ was performed by the authors of Reference (|130f ). who 
employed SU(2) L xSU(2)# %PT with explicit A isobar degrees of freedom and computed 
all contributions to 0(Q 3 ). These contributions arise from one-loop diagrams of the type 
illustrated in Figure [TT] Naively, one expects the loop contributions to be suppressed 
by powers of (Q/A x ) k for k = 2,3, where A x = 4iTF n ~ 1 GeV and Q is either m n or 
wa — m N . In the case of h\, however, the one-loop contributions receive logarithmic and 
fortuitous numerical enhancements, leading to the renormalized coupling 

hi = 0.5/i 1 + 0.25/4 - 0.24/ia + 0.08/^ , (43) 

where h\ is the bare nNN PV Yukawa coupling, h/\ is the analogous PV nNA Yukawa 
coupling, h\ parameterizes the isovector irNN PV derivative coupling 

C\ NN = v^N>f>ftN (n + D,n- - 7r~D^ + ) +■■■ (44) 

7T 

and parameterizes analogous PV 7riVA derivative couplings. 

The corrections appearing in Eq. (J43j) are only those having non-analytic quark 
mass dependence m q \nm q or m^ 2 (in the ira = wi n limit) and are uniquely identified 
with chiral loops 12 . As discussed in Reference (|l30l ). terms of this form cannot arise in 
quark model matrix elements of the hadronic weak Hamiltonian that were used in the 
analysis of DDH. Moreover, the "sum rule" contribution to h\ that DDH derived from 
AS = 1 decays using 811(6)^ symmetry relied on tree-level symmetry relations that do 
not contain m g -dependent symmetry breaking effects generated by chiral loops. Thus, 
it appears unlikely that DDH benchmark estimates for h\ fully reflect the impact of its 
quark mass dependence. Interestingly, the magnitudes of the coefficients of the various 
terms in Eq. (|43|) are comparable, allowing for possible cancellations between them that 
could reduce the magnitude of h\ from the DDH "best value" . In principle, a study of h\ 
on the lattice could allow one to identify the unknown constants appearing in Eq. (f4"3*|) 
by varying both m q and the number of colors. 

At present, carrying out such an analysis with unquenched QCD is prohibitively 
expensive. An alternative approach that allows one to extrapolate lattice computations 

11 The advent of chiral quarks has reduced the range over which an extrapolation must be performed. 
12 Terms that are analytic in m q can be absorbed into corresponding terms in the Lagrangian 
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Figure 11: One loop contributions to renormalized, PVirNN Yukawa coupling. 



to the domain of the physical, light quarks is to vary the valence and sea quark masses 
independently - a technique known as partial quenching. In order to identify the rng ence 
and m^ ea dependence analytically from one-loop computations, one mus t generalize stan- 
dard %PT to the corresponding partially-quenched effec tive theory (|l3lh . A computation 
using this framework has been performed in Reference (|l32l ). leading to an analogous ex- 
pression to that of Eq. (}4*3*|) that gives the non-analytic valence and sea quark mass 
dependence in the partially-quenched theory. An extrapolation to the chiral domain also 
requires inclusion of analytic terms that can be obtaine d fro m the effective Lagrangian 
and that were not written down explicitly in Reference (jl32h . 

As of this writing, the first-principles QCD analysis of h\ has not advanced beyond 
the analytic work of references HIS). Given the new experimental efforts, the 
time is clearly right for an investment in a lattice computation of h\. On a longer 
term horizon, one would also hope to see lattice predictions of the constants Cj and Ci 
that parameterize the 0(Q) short-distance, PV four-nucleon operators discussed above. 
Obtaining such computations will likely depend on progress in lattice calculations of the 
low-energy, strong NN interaction. A program aimed in this direction is being carried 
out by Savage and collaborators (M.J. Savage, private communication). 

In the absence of first principles QCD computations of the PV constants, one may 
look to nucleon model calculations for guidance as to their magnitudes. The expectations 
derived from the DDH SU(6) TO /quark model treatment (0) has been discussed above. 
Alternate approach es ha ve recently been used to predict h\. In the three-flavor Skyrme 
model of Reference (jl33h . the dominant contribution arises from terms in the AS = 0, PV 
four-quark Hamiltonian that contain strange quark bilinears such as s^^su^^^u. These 
terms sample the "kaon cloud" that arises from the Wess-Zumino action via rotations 
of the chiral soliton in SU(3) space. The resulting prediction is 2g w ^ h\ i$ 3.4<jv. 
The corresponding prediction in the two-flavor Skyrme model is considerably small er. A 
prediction for a somewhat larger value has been obtained using QCD sum rules (134), 
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wherein one expands nucleon correlators in a pion background in terms of various quark 
and gluon condensates. Values for the latter - such as (qir a '-f5q} n - are taken from other 
nucleon properties, such as the strong nNN coupling, leading to h\ ~ 8g n . Both sets of 
model predictions are roughly consistent with expectations based on NDA as well as the 
updated "best values" obtained from the SU(6)„,/quark model approach (jfil llih 



3.4 Experimental Prospects 

In the foregoing sections, we have emphasized the need for a complete set of precise 
PV measurements in the NN and few nucleon systems that can be cleanly interpreted 
in terms of constraints on PV LECs. As summarized in Table H3 we have currently 
two significant measurements in hand carried out with low energy proton beams (A p v 
and A p z a ), and two currently underway at existing neutron facilities (A ^ a nd d(j) na /dz). 
As noted earlier, there is also an existing pd asymmetry measurement (J991 ) which poses 
additional theoretical challenges accounting for the angular dependence in both elastic 
scattering and breakup channels, that remains to be analyzed in a common framework. 

Some immediate prospects for improvement in neutron beam measurements will 
take advantage of the superior features of a high intensity pulsed beam facility, the Spal- 
lation Neutron Source (SNS), currently under construction at Oak Ridge, Tennessee. The 
SNS is anticipated to provide the world's most intense beams of cold pulsed neut rons , 
approaching or even surpassing the time averaged intensities of cw reactor sources (jl35h 
by 2008. A crucial advantage of pulsed beams over reactor sources for precision PV ex- 
periments results from the introduction of new diagnostic capabilities via time-of-flight 
analysis of the neutron energy - especially important for reducing 7-ray backgrounds and 
systematic error diagnosis as well as neutron polarization diagnostics. In addition, the 
construction of a new facility allows for the incorporation of technological advances in 
neutron guide instrumentation, particularly the use of high efficiency bent supermirror 
transport guides which eliminate direct line of sight between the apparatus and the cold 
moderator, thereby reducing 7-ray and neutron backgrounds without significant loss of 
beam flux. The SNS is b uilding a new dedicated beamline for Fundamental Neutron 
Physics (FnPB) (jl35l . Il3fi| ) which is optimized to the needs of a suite of precision exper- 
iments in HWI and neutron beta decay. 



Neutron Capture Gamma Asymmetry Measurements 

As discussed earlier, the LANSCE phase of the ongoing NPDGamma experiment to 
measure A^ will be statistics limited at the 10~ 7 level, while the expected asymmetry 
based on the NDA and DDH 'best value' estimates for h\ is A^ = — 5 x 10 -8 . NPDGamma 
i s exp ected to be a key element of the initial SNS fundamental neutron physics program 
(137). At the time of writing, the apparatus has been commissioned, and systematic 
errors have been extensively studied at LANSCE, awaiting installation of the liquid 
parahydrogen target. A conclusion of these studies is that the apparatus is ready to 
make a measurement of A^ with a statistical error of 1 x 10~ 8 . Systematic errors are 
expected to be at or below the 10~ 9 level; however, additional running time to measure 
the potential false asymmetry from aluminum vacuum windows is needed to reach the 
ultimate experimental precision that can be achieved at the SNS. 
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The SNS FnPB beamline has been extensively studied via Monte Carlo simula- 
tions, and reasonably conservative flux estimates indicate that a measurement of at 
the 1 x 1CT 8 level should be possible in approximately 5000 hours of running on the 
new beamline once the SNS reaches 1.4 MW operation. The experiment can be moved 
with only minimal changes required to the apparatus, since the beam conditions (apart 
from increased flux) will be quite similar to those at LANSCE. It is planned to use the 
NPDGamma apparatus to commission the new FnPB beamline when it comes on line, 
beginning in 2008. In addition to the higher beam flux, two significant improvements 
in experimental conditions at the SNS are anticipated that could further improve the 
experimental precision. Gamma background reduction will be achieved in part 
suit of the SNS curved neutron guide, and in part with the incorporation of additional 
lead shielding upstream of the main detector array. Higher beam polarization should 
also be possible, via a combination of increased laser pumping power and/or the use of 
spectrally narrowed lasers for optical pumping of the 3 He spin filter cell, which have been 
demonstrated to produce up to 75% polarization in bench tests as compared to 4 0-55% 
3 He polarization routinely achieved during extended running in FP12 at LANSCE (|137l ). 

As noted earlier, the PV gamma asymmetry A* in the reaction n + d — > t + 7 
provides a complementary window on the AS = HWI to its counterpart in the np 
system. A Letter of Intent has been submitted to th e SN S to develop this experiment 
as a logical follow up to the NPDGamma experiment (|l38f ). As indicated in Table EH A^ 
displays a much larger sensitivity to the PV LEC's - including a ~ 30 times stronger 
sensitivity to h\ - than does A^. Based on the various estimates for the PV LECs 
(Table |2J), one could expect A^ to be one to two orders of magnitude larger than A^. 

In principle, the nd asymmetry measurement can be performed using most of the 
components of the NPDGamma apparatus, with an obvious exception of the target. The 
experiment is technically much more challenging than the np case due to the much smaller 
nd capture cross section; consequently, most of the neutrons will scatter out of the target 
rather than being captured to produce the gamma rays of interest. A room temperature 
liquid D 2 target is under consideration, with a length optimized between two compet- 
ing factors - the very small nd capture cross section, and the desire to avoid significant 
neutron depolarization in the target. A novel target vessel and shielding scheme will 
be required to absorb the scattered neutrons without producing significant background 
gamma rays. The interactions of polarized cold neutrons in D2O are at present not well 
understood, and a program of detailed simulations and test measurements will be carried 
out in order to optimize the design of the experiment. Other target possibilities include 
cold solid orthodeuterium or solid ortho-D 2 in an effort to minimize depolarization by 
slowing down the neutron beam, thereby increasing the capture probability (W. M. Snow, 
private communication) . An earlier measurement of A 1 was car ried out at ILL and found 
a result consistent with zero: A\ = (4.2±3.8) x 10" 6 (|Sll40h : the goal of the proposed 
SNS measurement is to reach a sensitivity 4 x 10~ 7 . Many of the systematic effects are 
similar to those for the np experiment, where they have been extensively studied. It 
should be noted that the absolute tolerance for systematic errors is relaxed for the nd 
experiment, since the absolute precision goal is more than an order of magnitude less 
stringent than that for the np experiment. 
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Neutron Spin Rotation Measurements 

As for the neutron capture asymmetry measurements, the higher neutron flux anticipated 
for the SNS as well as the pulsed nature of the beam offer compelling advantages for 
improving the precision of the na PV spin rotation measurement, as well as the possibility 
of carrying out experiments on np and perhaps even nd sp in rotation. To date, a Letter 



of Intent has been submitted for the na experiment (|14lh . anticipating a measurement 



accuracy of 1 x 1CT 7 rad/m in 12 months of data taking, which would represent a factor 
of three improvement over what is currently expected at NIST, with improved systematic 
error control. It is likely that the na spin rotation experiment will run at the SNS quite 
early in the FnPB experimental program. T he co llaboration intends in the longer term to 
pursue np spin rotation at the SNS as well (|l42l ). The sensitivity of the na spin rotation 
experiment to the five PV LEC's in the pionless EFT is given in table El ( predictions in 
the DDH meson-exchange model are available for the np and na cases, as mentioned in 



The basic spin rotation experimental technique has been outlined in Section 12.2.31 
recall the crucial role of the c 7r coil' in Figure which is used to suppress the effect 
of the much larger spin precessions due to residual magnetic fields inside the apparatus. 
Since previous experiments have been carried out at reactor facilities with polychromatic 
cold neutron beams, the best one could do was to optimize the field in the ir coil to rotate 
the neutron spins on average by 180° between the upstream and downstream segments 
of th e apparatus. Nearly a factor of two improvement can be obtained at a pulsed source 
(jl42f ) by ramping the ir coil current as a function of time of flight to provide a 180° 
spin rotation for all neutron velocities. In addition, one can take advantage of the fact 
that the PV spin rotation is independent of neutron velocity, whereas scattering and 
magnetic field effects are in general energy and thus velocity dependent for diagnosis 
of systematic effects. For the hydrogen and possible deuterium target measurements, 
neutron depolarization must be avoided by preparing cryogenic liquid targets in spin- 
selected molecular states - paramo lecular for hydrogen and orthomolecular for deuterium 
- an additional complication, particularly considering the requirement for constant cycling 
of the liquid between "front" and "rear" target locations. For neutron energies below 15 
meV, depolarization does not occur in parahydrogen. Less is known about the case for 
orthodeuterium, but a recent measurement at PSI found negligible depolarization of cold 
neutrons on orthodeuterium for a 4 cm target (W. M. Snow, private communication), 
which is encouraging news for the prospect of a deuterium spin rotation measurement. 

4 Beyond Hadronic Parity Violation 

Throughout this article, we have emphasized the significance of hadronic PV as a probe of 
the strangeness conserving HWI. Here, we comment briefly on its implications for other 
processes. We have already encountered one such process - PV electron scattering - 
discussed in Section 12.11 There, contributions from hadronic PV to the P V asymmetries 
constitutes a theoretical background that one must compute reliably in order to extract 
information on other quantities of interest, such as the strange quark form factors. 

A second illustration concerns neutrinoless double /3-decay (0i//5/3). This process has 




Section EH - 
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considerable current interest, as it provides the only known way to determine whether or 
not neutrinos are Majorana fermions (for recent reviews, see e.g., references ( 143l . 144j ).) 



In addition, one also hopes to use Ov/3/3 to determine the absolute scale of neutrino mass, 
complementing what we know about m u from tritium /3-decay and about neutrino mass 
differences from oscillation studies. 

Unfortunately, the latter use of QvfiP is complicated by possible contributions to the 
rate from the exchange of heavy Majorana particles, such as a heavy Majorana neutrino 
or the neutralinos of supersymmetry. In order to determine the absolute scale of m v , 
one must know the amplitude Ah for heavy particle-exchange contributions, which can 
be comparable in magnitude to the amplitude At, fo r light Majorana neutrino-exchange. 
From simple dimensional arguments, one has (145) 



^l^MML , (45 ) 

A L A 5 m^ 

where mpp is the effective mass of the light Majorana neutrino having typical virtuality 
P ~ (50 MeV) 2 and A is the mass scale associated with the heavy Majorana particles. 
Given what we know about AmJ and the neutrino mixing matrix elements that help 
determine mpp, the ratio in Eq. ()45|) can be 0(1) for A of order 1 TeV. Thus, it is 
important to analyze the possible heavy particle contributions with theoretical clarity. 

Rece ntly, an EFT approach for doing so was developed by the authors of Refer- 
ence (Il46f ). At the lepton-quark level, the effective operators for heavy, Majorana particle 
exchange factorize into products of four-quark and two-lepton operators. The four-quark 
operators are analogous to those entering the AS = HWI, differing only in their re- 
spective representations in chiral SU(2). Consequently, the mapping of these operators 
onto effective, hadronic operators involving nucleon and pion degrees of freedom is simi- 
lar to the one used in obtaining the EFT for hadronic PV. Thus, the study of hadronic 
PV should provide insights into the EFT for heavy particle contributions in Qvf3f3. We 
consider two aspects of this correspondence in particular: 

(i) In contrast to the situation for hadronic PV, there does not exist a program of 
few-body experiments from which one can determine the operator coefficients for 
the Oz//3/3 EFT. These coefficients must be computed theoretically. Comparisons 
of analogous computations of the LECs for hadronic PV with experimental values 
should provide new guidance for obtaining reliable computations in the Oz//3/3 case. 
Indeed, as argued above, hadronic PV provides a unique tool for learning how the 
strong interaction dresses four-quark weak interactions into hadronic operators and 
amplitudes. 

(ii) Again in contrast to hadronic PV, the only systems in which 0uf3f3 can occur 
are heavy, complex nuclei. While there has been considerable recent progress in 
understanding how an EFT power-counting of operators translates into a power- 
counting of few-body matrix elements, the situation with complex nuclei is less 
clear. In order for the EFT to provide realistic guidance to the size of heavy 
particle contributions to 0z//3/3 transition matrix elements, we need to know how 
well operator power counting works for nuclear matrix elements. 
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In principle, hadronic PV can provide useful insights into this issue. Since the 
lowest order PV interaction can be determined from a program of few-body exper- 
iments as outlined above, this interaction can then be used as a known probe of 
nuclear PV observables, such as the nuclear anapole moment or the PV 7-decays 
of p-shell nuclei. To the extent that the lowest order PV interaction suffices to 
yield successful descriptions of these nuclear PV observables, we would have evi- 
dence for the applicability of the EFT to nuclei. To the extent that it does not, we 
would conclude that many-body renormalization of the lowest-order weak effective 
interaction is substantial and that high-momentum components of nuclear wave- 
functions play a more important role than one might naively expect. Either way, 
the implications for 0vf3f3 would be important. 

5 Summary and Conclusions 

The quest to explain the manifestations of weak interactions between quarks in strongly- 
interacting systems remains an important piece of "unfinished business" for Standard 
Model physics. In both the AS = 1 and AS = sectors , the non-perturbative character 
of low-energy QCD has been the stumbling block. The AS = 1 decays of hyperons 
in particular seem to elude explanation using the standard symmetry and effective field 
theory approaches that have been so successful in treating low-energy strong interactions. 
Whether these puzzles simply reflect the active participation of the strange quark with 
its mass of order the QCD scale, or some other dynamics peculiar to hadronic weak 
interactions, is unknown. Our hope is that through experimental studies of the AS = 
HWI with PV observables and through their theoretical interpretation within the EFT 
framework that makes for the closest possible contact with QCD, we will gain new insights 
into the low-energy weak interactions of the three lightest quarks. In this review, we hope 
to have sketched a useful roadmap for future progress in this direction. 

As we have emphasized throughout this article, the forefront in this endeavor lies in 
the arena of single-nucleon and few-body nuclear systems. Experimental developments 
have paved the way for completion of precise measurements of the O(10~ 7 ) PV effects 
in such systems, and plans are underway for several new few-body experiments. The 
immediate theoretical challenge is to compute the few-body PV observables using the 
EFT framework, allowing one to extract robust values for the low energy constants from 
experiment. In the longer term, we would like to determine the extent to which these 
results are consistent with Standard Model-based expectations, both from the standpoint 
of symmetries and from first principles lattice calculations. In light of the substantial ex- 
perimental and theoretical progress in the field as well as the significant new opportunities 
made possible by this progress, we are optimistic that this effort can be successful. 
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